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ABSTRACT 

We present radial velocities and chemical abundances for nine stars in the old, distant open clusters Be 18, 
Be 21, Be 22, Be 32, and PWM 4. For Be 18 and PWM 4, these are the first chemical abundance measure- 
ments. Combining our data with literature results produces a compilation of some 68 chemical abundance 
measurements in 49 unique clusters. For this combined sample, we study the chemical abundances of open 
clusters as a function of distance, age, and metallicity. We confirm that the metallicity gradient in the outer disk 
is flatter than the gradient in the vicinity of the solar neighborhood. We also confirm that the open clusters in 
the outer disk are metal-poor with enhancements in the ratios [a/Fe] and perhaps [Eu/Fe]. All elements show 
negligible or small trends between [X/Fe] and distance (< 0.02 dex/kpc), but for some elements, there is a hint 
that the local {Rgc < 13 kpc) and distant (Rgc > 13 kpc) samples may have different trends with distance. 
There is no evidence for significant abundance trends versus age (< 0.04 dex Gyr" 1 ). We measure the linear 
relation between [X/Fe] and metallicity, [Fe/H], and find that the scatter about the mean trend is comparable to 
the measurement uncertainties. Comparison with solar neighborhood field giants shows that the open clusters 
share similar abundance ratios [X/Fe] at a given metallicity. While the flattening of the metallicity gradient and 
enhanced [a/Fe] ratios in the outer disk suggest a different chemical enrichment history to the solar neighbor- 
hood, we echo the sentiments expressed by Friel et al. that definitive conclusions await homogeneous analyses 
of larger samples of stars in larger numbers of clusters. Arguably, our understanding of the evolution of the 
outer disk from open clusters is currently limited by systematic abundance differences between various studies. 
Subject headings: Galaxy: abundances, Galaxy: disk, Galaxy: open clusters and associations: general 



1. INTRODUCTION 

Ou r Galaxy's open clusters are valuable tools to study the 
disk ((Friel 1995). Accurate homogeneous distances and ages 
can be measured for samples of open clu sters that span a wide 
range in parameters (Sala ris et al.ll20 04). Additionally, metal- 
licity estimates of open clusters can be readily obtained from 
a variety of methods (albeit with differing degrees of accu- 
racy) thereby enabling studies of the structure, kinematics, 
and chemistry of the disk as well as any temporal variations 
of these properties. 

The atmospheres of low-mass stars retain, to a great extent, 
the chemical composition of the interstellar medium at the 
time and place of their birth. The nucleosynthetic yields of the 
chemical elements depend upon stellar mass and metallicity. 
Therefore, measurements of metallicity, [Fe/H], and chemical 
abundance ratios, [X/Fe], in stars in open clusters offer power- 
ful ins ight into the formation and ev olution of the disk (I Jane si 
1979t iFreeman & Bland-Hawthornl [20021: iFriel et all 120021: 
Bland-Hawthorn et al.ll2010l:lKobayashi & Nakasatoll201 lh . 

In recent times there has been considerable effort to under- 



The data presented herein were obtained at the W.M. Keck Observa- 
tory, which is operated as a scientific partnership among the California In- 
stitute of Technology, the University of California and the National Aero- 
nautics and Space Administration. The Observatory was made possible by 
the generous financial support of the W. M. Keck Foundation. 



stand th e evolution of th e outer G alactic disk (e.g.jHou et al. 
2000t IChiappini et all 1200 U l Andrievskv et all 12002c 
Daflon & Cunhal2004tlCosta et al.ll2004tlCescutti et alJl2007 



Magrini et alJl2009h . It appears that the metallicity gradient 
(i.e., [Fe/H] versus Galactocentric distance, Rgc) in the 
outer disk {Rgc > 13 kpc) is flatter t han the metallicity 



gradient in the so l ar neighborhood (e.g.. iTwarog et al.l 1 1 997 



Luck et al 



Yong et al 



2003 



2005 



Carraro et all 12004 ICarnev et all 12005b 
Bragaglia et all 120081 ISestito et al.l 12008 



Jacobs on et al T l2009HFriel et al.l 1201 Oh . Qualitatively similar 



behavior has now been found in external galaxies (e.g. . 
Worthev et al.ll200H iBresolin et al J 120091: iVlaiic et alJl2009l 



2011b suggesting that the processes which govern chemical 



evolution in the outer Galactic disk have also operated in the 
outskirts of other disk galaxies. 

Another intriguing result is evidence for enhanced 
ra/Fe ] ratios in th e oute r Galactic disk fro m open clus- 
t ers dCarraro et al.l 12004 lYong et ail 120051) . field stars 
Carney et alj l2005bt iBensbv et alj 1201 lh , and Cepheids 



Yong et al.ll2006h~ Such abundance ratios in objects spanning 



a large range in ages can be explained by infall of pristine gas 
in the outer disk and/or vigorous star formation. Either expla- 
nation requires that the outer disk has experienced a signifi- 
cantly different star formation history compared to the solar 
neighborhood and this has important implications for the for- 
mation and evolution of the Galactic disk. 
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Be 18 (Kaluzny 1997) 
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FIG . 1. — Be 18 color-magnitude diagram using the data from IKaluznvl 
<1997t ). We distinguish our program stars by red squares. 
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FIG. 2. — Be 18 c olor-magnitude diagram using the data from 2MASS 
Skrutskie et al. 12006). We distinguish our program stars by red squares. 



However, not a ll studies of the oute r disk find enhanced 
[a/Fe] ratios (e.g jBragaglia et ai]|2008i) . While one explana- 
tion for the discrepancy may be systematic differences in the 
analyses, studies of additional open clusters in the outer disk 
are necessary to clarify the situation. Therefore, to further ex- 
plore the formation and evolution of the outer Galactic disk, 
we analyze the chemical abundances of a new sample of open 
clusters. This is the fourth, and final, paper in our series on 
element al abundance ratio s in stars in the outer Galactic disk. 
Paper I dYong et alj 120051) co ncentrated upon open clusters, 
Paper II dCarnev et alj 2005bt) was d edicated to field red gi- 
ants, and Paper III dYong et aUl2006h focused upon Cepheids. 
In this paper, we present radial velocities, metallicities [Fe/H], 
and element abundance ratios [X/Fe] for the open clusters Be 
18, Be 21, Be 22, Be 32, and PWM 4. 

2. PROGRAM STARS AND OBSERVATIONS 

2.1. Target Selection 

We searched the WEBDA databas^3 and the recent litera- 
ture to identify additional old (> 2 Gyr), distant (Rgc > 10 
kpc) open clusters. The following open clusters were identi- 
fied as suitable candidates: Be 18, Be 21, Be 22, Be 32, and 
PWM 4. Following our approach in Paper I, individual stars 
were selected from optical and infra-red color-magnitude di- 
agr ams (see Figures [D to [TOl and Tables [T] and |2). For Be 
32. lD'Orazi et alJd2006l) measured radial velocities for a large 
sample of stars and thus we were able to select likely members 
for this cluster. 
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FIG. 3. — Be 21 color-magnitude diagram using the data from ITosi et alj 
11998ft . We distinguish our program stars by red squares. 



2.2. Observations and Data Reduction 

We used the HIRES dVogt etalJll994h spectrograph on the 
Keck I Telescope on 12 November 2006 and 27 December 
2006. On both occasions we used the red cross-disperser and 
the B5 decker which has a slit width of 0."86 and a slit length 

2 http://www.univie.ac.at/webda/ 



of 3. "5. On the November run, we binned the CCD 3x1 
(spatial x spectral) and on the December run we binned the 
CCD 2x1. The observing routine included 20 quartz lamp 
exposures for flat fielding as well as 20 zero second exposures 
for "bias" frames. ThAr frames were taken during each night 
(three on 1 2 November and eight on 27 December) and two 
radial velocity standards were observed on each night. The 
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TABLE 1 
Observed Clusters 



Cluster 


R.A. a (J2000) 


Decl. a (J2000) 


£» 


V 


[Fe/Hf 


SV b 


MAI C 


d d 


d e 


Berkeley 18 


05:22:12 


+45:24:00 


163.63 


+5.02 


+0.02 


2.3 


5.69 


5.8 


5.4 


Berkeley 21 


05:51:42 


+21:47:00 


186.84 


-2.51 


-0.83 


1.6 


2.18 


5.0 


6.2 


Berkeley 22 


05:58:24 


+07:50:00 


199.88 


-8.08 


-0.30 


2.1 


4.26 


7.7 


6.2 


Berkeley 32 


06:58:06 


+06:26:00 


207.95 


+4.40 


-0.58 


2.4 


5.91 


3.1 


3.2 


PWM4 


23:50:55 


+62:19:15 


115.96 


+0.27 




3.0 


7.00 f 


7.9 f 


7.2 



a Taken from WEBDA. 

Taken fronrSalaris et 
c Taken from Salaris et al. i 2004). in Gyrs. 
d Distance in kpc from lSalaris et aljt20t3) . 
e Distance estimate obtained using red clump stars (this paper). 
^ Taken from Ortolani et al. ( 2005) 



TABLE 2 
Photometric Data 



Star R.A. a (J2000) Decl. a (J2000) V B-V V-I c J-K" Reference 



Bel8 532 


05 22 13.7 


+45 25 46.8 


15.99 


1.48 


1.68 


12.226 


0.877 


1 


Bel8 1006 


05 22 16.2 


+45 27 28.5 


16.03 


1.47 


1.70 


12.224 


0.885 


1 


Bel8 1163 


05 22 13.8 


+45 27 58.9 


15.79 


1.63 


1.86 


11.561 


1.010 


1 


Bel8 1383 


05 22 10.7 


+45 28 49.4 


15.07 


1.75 


1.94 


10.631 


1.045 


1 


Be21 T50 


05 51 42.3 


+21 48 45.4 


15.87 


1.71 


2.07 


11.142 


1.062 


2 


Be21 T51 


05 51 42.0 


+21 48 02.8 


15.69 


1.70 


2.06 


10.963 


1.095 


2 


Be22 414 


05 58 25.9 


+07 46 11.4 


15.76 


1.87 


2.02 


11.020 


1.157 


3 


Be22 643 


05 58 26.9 


+07 45 26.1 


14.31 


2.29 


2.65 


8.369 


1.405 


3 


Be32 16 


06 58 06.9 


+06 25 56.5 


13.61 


1.06 




10.970 


0.660 


4 


Be32 18 


06 58 13.8 


+06 27 54.9 


13.71 


1.06 




11.069 


0.640 


4 


PWM4 RGB1 


23 50 57.4 


+62 20 03.2 


16.19 


2.21 




9.876 


1.446 


5 


PWM4 RGB 2 


23 50 58.3 


+62 19 18.0 


16.91 


1.99 




11.185 


1.305 


5 


PWM4RC1 


23 50 53.9 


+62 19 07.7 








13.015 


1.142 




PWM4 RC2 


23 50 53.2 


+62 19 33.8 








13.087 


1.087 





REF ERE NCES. — 1 =IKaluznvl {19971) : 2 = ITosi et alj fi99"8t) : 3 = IKaluznvl fl994h : 4 = IKaluznv & MazuJ 
1 1991); 5 = Ortolani et al. <2()05|) 



11 2MASS (Skrutskie et al. 2006:) coordinates and photometry. 



spectroscopic data were reduced using the IRAF0 packages as 
described in Paper I. From the reduced ThAr frames, we mea- 
sured a spectral resolution of R = 47,000 (cr = 3000) on each 
of the three detectors, and this value is close to the expected 
value of R = 48,000. Since we did not bin the CCD in the spec- 
tral dimension, the data are oversampled with approximately 
4.8 pixels per resolution element. The wavelength coverage 
was from 4000A to 8500A. There are gaps in the spectra be- 
tween the blue and green CCDs (~ 5400 A to 5440 A) and 
between the green and red CCDs (~ 7000A to 7075A). Fur- 
thermore, the free spectral range exceeds the CCD coverage 
beyond - 6400A. 

Exposure times and resulting signal-to-noise ratios (S/N) 
for the program stars are given in Table [3] (Note that the S/N 
is given per pixel at 6050A and that we have 4.8 pixels per 
resolution element. To obtain S/N ratios per resolution ele- 
ment, multiply our values by 2.19.) Our observing program 
consisted of relatively short exposures of some stars in order 
to measure their radial velocities as well as (multiple) longer 
exposures of other stars to obtain high S/N ratios for detailed 

3 IRAF (Image Reduction and Analysis Facility) is distributed by the Na- 
tional Optical Astronomy Observatory, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. 



chemical abundance analysis. For the longer exposures, pre- 
liminary reduction was performed in real-time in order to es- 
timate the radial velocities. Shorter exposures were taken of 
four candidate red clump giants (two stars in each of Be 18 
and PWM4) and one red giant in PWM4, and the radial veloc- 
ities were measured during the analysis. Examples of reduced 
spectra are shown in Figure [TT] 

2.3. Radial Velocities 

Radial velocities were measured by cross-correlating the 
program stars with radial velocity standards. For 12 Novem- 
ber 2006, our radial velocity standard was HD 90861 with 
a heliocentric radial velocity of 36.3 km s" 1 (this value was 
used in Paper I and was taken from the Astronomical Almanac 
1999). For 27 December 2006, our radial velocity standard 
was HD 1 10885 with a heliocentric radial velocity of -47.89 
km s" 1 dCarnev et alj|2003l) . 

The cross-correlation typically involved some 30 orders per 
star and the dispersion about the mean velocity ranged from 
0.1 to 0.7 km s . However, we found that on 12 November 
2006, the velocity "jumps" by ~ 2 km s" 1 within the consec- 
utive series of spectra taken for Be21-T51. On closer inspec- 
tion of the data, we found that one of the four spectra obtained 
has an anomalous radial velocity. The wavelength solution 
for this spectrum comes from a different ThAr frame than for 



4 



10 - 



12 



14 



1 i 1 1 1 i 1 1 1 + 1 

Be 21 (r < 150") 



-i — i — i — i — r- 



+ 
t 



+ 

+ + + 





+ 


+ 

+ 


+ 




++ 


+ 








+ + + 








+ + 




+ 


+ 




++ + 
+ $+ 


++ + 


++ 


+ 




% +* 
+ ++ 


++ + 


+ 







_l I I I I I I I I I I I I I I I I I I I I I L_ 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 
J-K (2MASS) 

FIG. 4. — Be 21 c olor-magnitude diagram using the data from 2MASS 
Skrutskie et al. (2006). We distinguish our program stars by red squares. 
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FIG. 6. — Be 22 c olor-magnitude diagram using the data from 2MASS 
Skrutskie et al. (2006). We distinguish our program stars by red squares. 
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FIG . 5. — Be 22 color-magnitude diagram using the data from IKaluznvl 
(1994!). We distinguish our program stars by red squares. 



FIG. 7. — Be 32 color-magnitude diagram using the data from 
Kaluzny & Mazur 1 1991). We distinguish our program stars by red squares. 



the other three spectra. On this observing run, we only took 
three ThAr frames, one each at the start and end of the night 
and only one in the middle of the night. At face value, it 
would appear that there has been a shift during the course of 
the night and that the radial velocities on this night probably 
have a systematic uncertainty of ~ 2 km s" 1 . This issue was 
noted shortly after the November run and more frequent ThAr 
frames, eight in total, were taken on the December run. For 
Be21-T51, measurement of the telluric absorption line wave- 
lengths suggests that the three spectra with consistent radial 



velocities give the correct value. Another point worth noting 
is that Bel8-1163 was observed on both runs. For this star, 
we find that the radial velocity from 12 November 2006 is 
~1 km s" 1 higher than the values obtained from 27 December 
2006. This could be explai ned by either radial velocity varia- 
tion (and/or jitter - e.g., see lCarneV et al. 2008|and references 
therein) or systematic errors in the radial velocity standards. 
The bottom line is that our radial velocities may be uncertain 
by ~1 km s" 1 between the two observing runs. 
For Be 18 and PWM4, there are no previous radial velocity 



5 



8 



10 - 



12 



14 - 



Be 32 (r < 180") 



+ 

+ + 

+ 

+ 

+ + 



+ 
+ 



+ 

+ + + + , 



_i I i i i I i i i I i i i I i i i I i i i I i_ 

0.0 0.2 0.4 0.6 0.8 1.0 
J-K (2MASS) 

FIG. 8. — Be 32 c olor-magnitude diagram using the data from 2MASS 
ISkrutsk ie et al. (2006j). We distinguish our program stars by red squares. 
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FIG. 10. — PWM 4 color-magnitude diagram using the data from 2MASS 
Skrutskie et al. (2006). We distinguish our program stars by red squares. 
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PWM 4 (Ortolani et al. 2005) 




TABLE 3 
Spectroscopic Observations 



2.5 



FIG. 9. — PWM 4 color-magnitude diagram using the data from 
I Ortolani et al. (2005). We distinguish our program stars by red squares (the 
two red clump stars are not included in their photometry). 



studies to our knowledge. As noted, Bel8-1 163 was observed 
on both runs and there was a 1 .0 km s" 1 radial velocity offset 
between the two runs. Without knowing which of the two val- 
ues may be spurious, we include all values for all stars when 
computing the average cluster radial velocity for Be 18 of -5.5 
±1.1 km s" 1 (cr = 2.2 km s" 1 ). For PWM4, we find an aver- 
age cluster radial velocity of -125.9 ± l.Okms" 1 (cr = 2.0km 
s" 1 ). For both clusters we studied four stars. 
For Be21 T50, our average radial velocity is -1.4 km s" 1 



Star 


Exposure Time (s) 


HJD - 2,450,000 


S/N a 


v rad 


Bel8-532 




600 


4052.9209 


22 


-3.6 


Be 18- 1006 




500 


4052.9325 


20 


-8.3 


Bel8-1163 




500 


4052.9414 


24 


-3.1 


Bel8-1163 




2400 


4097.7714 


48 


-4.1 


Bel8-1163 




2401 


4097.7998 


51 


-4.4 


Bel8-1163 




2400 


4097.8282 


50 


-4.4 


Bel8-1383 




1200 


4052.9535 


48 


-6.5 


Bel8-1383 




1800 


4052.9717 


58 


-6.1 


Bel8-1383 




1500 


4052.9916 


54 


-6.1 


Be21-T50 




2400 


4097.9890 


53 


-1.3 


Be21-T50 




2401 


4098.0174 


54 


-1.6 


Be21-T50 




2401 


4098.0458 


51 


-1.2 


Be21-T51 




2064 


4053.0212 


53 


+0.9 


Be21-T51 




2400 


4053.0535 


57 


-1.1 


Be21-T51 




2400 


4053.0818 


58 


-0.9 


Be21-T51 




600 


4053.0998 


28 


-0.9 


Be22-414 




2400 


4097.8690 


51 


+93.1 


Be22-414 




983 


4097.8893 


31 


+93.0 


Be22-414 




2400 


4097.9184 


51 


+93.4 


Be22-414 




2533 


4097.9487 


52 


+93.7 


Be22-643 




1200 


4053.1184 


72 


+90.4 


Be22-643 




1200 


4053.1328 


69 


+90.4 


Be32-16 




600 


4098.0674 


61 


+105.9 


Be32-16 




600 


4098.0751 


61 


+105.9 


Be32-16 




600 


4098.0827 


62 


+105.8 


Be32-18 




600 


4098.0910 


58 


+105.1 


Be32-18 




600 


4098.0985 


54 


+105.1 


Be32-18 




901 


4098.1080 


64 


+105.1 


PWM4-RGB 1 




2002 


4052.7283 


42 


-125.4 


PWM4-RGB 1 




2700 


4052.7631 


48 


-124.8 


PWM4-RGB 1 




2700 


4052.7953 


51 


-124.8 


PWM4-RGB 1 




2700 


4052.8276 


51 


-124.6 


PWM4-RGB2 




600 


4052.8527 


17 


-128.8 


PWM4-RC1 




1800 


4052.8709 


17 


-125.4 


PWM4-RC2 




1800 


4052.8963 


16 


-124.3 


a Signal-to-noise ratio per pixel near 6050A. 

(""internal = 0-2 km s" 1 ) and this value agrees well with the aver- 
age cluster value of -0.6 km s -1 (cr = 2.9 km s -1 ) from Paper I. 
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FIG. 1 1 . — Sample region of spectra of the program stars. Lines used in the 
analysis are marked. 

For Be21 T51, our average radial velocity from all four spec- 
tra is -0.5 km s -1 (<7i nterna i = 0.9 km s" 1 ). Had we excluded the 
apparently anomalous +0.9 km s" 1 value, we would have ob- 
tained an average value of -1.0 km s" 1 (<7j nte rnai = 0.1 km s -1 ). 
Either value is in good agreement with the value measured for 
this star in Paper I, -1.0 ±1.0 km s . However, as noted in 
Paper I, our radial velocities for this cluster differ from those 
of iHill&Pasquinild 19991) who obtained +12.35 ± 0.6 km s _1 
based on four stars. 

For Be22, we find an average cluster radial velocity of 
+91.9 ± 1.5 km s" 1 (^internal =2.1 km s" 1 ) based on two 
stars. This value is in fair agreement with the value from 
iVillanova et all (120051) of +95.3 ± 2.0 km s" 1 (a = 2.8 km 
s" 1 ). 

For Be32 18, our average radial veloci ty of + 1 05 . 1 km s 1 is 
in good agreement with the values from lD'Orazi et all (12006b 
and |SestitoetiD (|2006j) of + 1 08 .07 km s" 1 (a = 4.12 km s _1 ) 
and +105.5 km s" 1 , respectively. We find an average cluster 
radial velocity of + 105.5 ± 0.4 km s" 1 ((Tintemai = 0.5 km s -1 ) 
based on two stars. Our value is in good agreement with pre- 
vious studies of th is cluster: +101 ±3 km s -1 ( a = 10 km s" 1 ) 
(IScott etal.1 119951) : +106.7 (a = 8.5 km s~ l ) dD'Orazi et al.l 
I2006J); +10 6.0 ± 1.4 km s" 1 (ISestito et al.ll2006h : + 105.4~± 
0.4 km s" 1 dFriel et al.ll20"Tol) . 

2.4. Distance Estimations 

We determined distances and interstellar reddening for the 
clusters following the prescription outlined in Paper I. From 



the 2MASS color-ma gnitude diagrams, we estimated the loca- 
tion of the red clump. ICarnev et ail (12005 al) derived a relation 
between the mean J-K color of the red clump and metallicity. 
Using our spectroscopic measurement for metallicity, this re- 
lation is used to derive the reddening E(J-K) [= Q.52E(B-V) 
dRieke & Le bofskv 1985)]. The luminosity of the red clump 
is Mk — -1.61 mag and is insens itive to metallicity in the 
regime [Fe/H] > -1.0 dAlvesll2000l) . We present our distances 
and reddening in Table |4] Our distance estimates are, in gen- 
eral, in good agreement with the WEBDA literature values. 
Such agreement is encouraging given that in some cases, iden- 
tifying the clump from 2MASS CMDs can be difficult. 

3. ELEMENTAL ABUNDANCES 
3.1. Stellar Parameters 

Stellar parameters are required to conduct a chemical abun- 
dance analysis and we obtained the values in the following 
manner. Following Paper I, initial estimates of T e ff come 
from the infra-red flux method metallic ity-dependen t color- 
tem perature relations for giant s tars by lAlonso et al.l (1 1999b 
and IRamfrez & Melendezl d2005l) . With distances, tempera- 
tures, and bolometric corrections, we can calculate logg val- 
ues if we assume a mass for the red giants. As in Paper I, we 
adopt a mass of 1M© and note that the surface gravity is not 
especially sensitive to the assumed value. 

These initial estimates for T e s and \ogg were then used 
as starting values for the spectroscopic stellar parameters. 
Equivalent widths (EW) were measured for a selection of Fe I 
and Fe II lines using routines in IRAF. We used the same line 
list as in Paper I. Mod el atmosphere s were computed using 
the ATLAS9 program dKuruczll 1993b. Using the LTE stellar 
line analysis program MOOG dSnedenlfl973l) . we computed 
abundances for each line based on the measured EW. (As in 
Papers I, II, and III, we used the 2002 version of MOOG.) 
We set the effective temperature, T e ff, by requiring that the 
abundances from Fe I lines showed no trend with lower ex- 
citation potential, i.e., excitation equilibrium. The microtur- 
bulent velocity, was determined by ensuring that the abun- 
dance from Fe I lines showed no trend with reduced equivalent 
width, log(WA/A). We adjusted the surface gravity, \ogg, un- 
til the abundance from Fe I lines agreed with the abundance 
from Fe II lines, i.e., ionization equilibrium. We also required 
that the derived metallicity, [Fe/H], was within 0.1 dex of the 
metallicity of the model, [m/H], otherwise the model atmo- 
sphere was re-computed with the appropriate metallicity. The 
initial photometric stellar parameters and the final spectro- 
scopic stellar parameters are given in Table [5] 

We note that there is good agreement between the photo- 
metric and spectroscopic values. Given the large reddening 
for the program clusters, the infra-red flux method tempera- 
tures may be regarded as being less reliable than spectroscopic 
values. Furthermore, while uncertainties in the distances esti- 
mates affect the photometric temperatures, the spectroscopic 
values are not affecte d. For PWM4-RGB1, t he pho tometric 
parameters, using the IRamfrez &~M elendezl d2005l) calibra- 
tion, differ from our spectroscopic values. Such a discrepancy 
may be attributed to the uncertainty in the reddening (we ob- 
tain E(B-V) = 0.99 compared to the WEBDA value of 0.62). 
With the exception of PWM4, that the two sets of temper- 
atures are in good agreement, in general, probably suggests 
that the reddening and distance estimates are reasonable. 

In practice, we explored discrete values for T e ff (every 25 K, 
i.e., 4525K, 4550K, etc) and logg (every 0.05 dex, i.e., 1.05 
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TABLE 4 
Red Clump Data 



Cluster 


<K> 


<J-K> 


[Fe/H] 


(J-K) 


E(J-K) 


E(B-V) 


(m-M) 


d (kpc) 


Rgc 


E(B-Vf 


d a (kpc) 


Berkeley 18 


12.21 


0.87 


-0.4 


0.53 


0.34 


0.66 


13.60 


5.2 


13.10 


0.46 


5.8 


Berkeley 21 


12.62 


0.93 


-0.3 


0.55 


0.39 


0.74 


13.98 


6.2 


14.22 


0.76 


5.0 


Berkeley 22 


12.57 


0.88 


-0.4 


0.53 


0.35 


0.68 


13.95 


6.2 


13.92 


0.70 


7.7 


Berkeley 32 


11.02 


0.65 


-0.3 


0.55 


0.11 


0.20 


12.56 


3.3 


10.97 


0.16 


3.1 


PWM4 


13.02 


1.06 


-0.3 


0.55 


0.52 


0.99 


14.29 


7.2 


12.91 


0.62 


7.9 



a Taken from WEBDA. 



TABLE 5 
Atmospheric Parameters 



Star 


7dT 




Tea 


l°g£ 


Tea 


l°g£ 


[m/H] 


6 


[Fe/H] 




Photometric" 


Photometric 




Spectroscopic 




Bel8-1163 


4666 


2.1 


4655 


2.1 


4500 


2.2 


-0.5 


1.20 


-0.46 


Bel8 1383 


4489 


1.7 


4495 


1.7 


4400 


1.9 


-0.5 


1.27 


-0.41 


Be21 T50 


4464 


1.7 


4468 


1.7 


4625 


1.9 


-0.3 


1.32 


-0.26 


Be21 T51 


4481 


1.7 


4452 


1.6 


4500 


1.7 


-0.3 


1.29 


-0.35 


Be22 414 


4340 


1.7 


4297 


1.7 


4350 


1.7 


-0.5 


1.18 


-0.41 


Be22 643 


3791 


0.6 


3752 


0.5 


3850 


0.2 


-0.5 


1.55 


-0.49 


Be32 16 


4977 


2.4 


5000 


2.4 


4875 


2.4 


-0.3 


1.07 


-0.39 


Be32 18 


4962 


2.4 


5007 


2.4 


4950 


2.7 


-0.3 


1.46 


-0.37 


PWM4 RGB 1 


4368 


1.4 


3942 


1.0 


3950 


0.5 


-0.3 


1.28 


-0.33 



a Estimates obtained using (i) the photo metric data in T abled (ii) the distance, reddening, 
and metallicity from Table|4] and (iii) the Ramirez & Melendez 12005) IRFM calibration. 
b Estimates obtained using (i) the photo metric data in Table [2] (ii) the distance, reddening, 
and metallicity from Table|4] and (iii) the Alonso et al. 1 1999) IRFM calibration. 
c Quantities are derived using the spectroscopic methods described in the text. 



dex, 1.10 dex, etc). We assumed that excitation equilibrium 
was satisfied when the slope between log e(Fe I) and lower ex- 
citation potential was < 0.004 and that ionization equilibrium 
was satisfied when | log e(Fe I) - log e(Fe II) | < 0.05 dex. The 
microturbulent velocity was considered satisfactory when the 
slope between log e(Fe I) and reduced equivalent width was 
< 0.004. We iterated until the criteria were simultaneously 
satisfied. During this process, we estimate that the internal 
uncertainties are T e ff ± 100K, \ogg ±0.3 dex, and £, ± 0.2 
km s" 1 . 

3.2. Elemental Abundance Analysis 

In Paper I, we relied upon spectrum synthesis to measure 
chemical abundances for individual elements. In that paper, 
the spectral resolution was R = 28,000 and we were con- 
cerned that blends might affect the abundances derived from 
an equivalent width analysis. In this paper, the spectral reso- 
lution R = 47,000 was sufficiently high that abundance mea- 
surements based on EW measurements from our high S/N 
spectra were regarded to be reliable. Therefore, we measured 
chemical abundances from spectrum synthesis as well as from 
equivalent width analysis, when possible, using MOOG in 
both cases. For both approaches, we took into account the 
effects of hyperfine splitting and/or isotopic splitting when 
necessary (e.g., Mn, Co, Rb, Ba, La, and Eu) using the same 
approach as in Paper I. The line list and solar abundances were 
the same as in Paper I. (See Table [6] for the atomic data and 
EW measurements.) In Tables 171 to [TT1 we present the chemi- 
cal abundances, [X/Fe], for the program stars where these val- 
ues represent the average abundances from EW analysis and 
from spectrum synthesis. We note that the two approaches 



gave quite similar results with an average difference in abun- 
dance ratio A[X/Fe] (EW - Synth) of 0.06 ± 0.01 dex (a = 
0.14). We did not find any significant global differences for 
a given line or a given species between the two approaches. 
(The [X/Fe] ratios are based on the individual stellar [Fe/H] 
rather than the cluster mean [Fe/H].) In Table [T2l we show 
the abundance dependences upon the model parameters for 
a representative star, Be21 T51, assuming that the errors are 
symmetric for positive and negative changes. 

3.3. Comparison With Literature 

iFriel et all (120101) conducted an examination o f the various 
system atic differences between th e work by (a) ISestito et al.l 
(120081) and lBragagliaet al.l (120081) (SB using their notation), 
(b) preliminary results from Th is Study (CY usin g their no- 
tation^!), an d ( c ) re-analysis by IFriel et all (120101) using the 
spectra employed by SB and CY. We refer the reader to 
their careful analysis that takes into account solar abundances, 
EW measurements, lo g;?/' values and atmospheric parame- 
ters. IFriel et all (120101) identify systematic abundance differ- 
ences between the various studies, but note that "it is difficult 
to identify the component sources of systematic differences". 

For completeness, we also include a comparison of mea- 
sured equ ivalent widths for Be3 2-18. For 67 lines in com- 
mon with B ragaglia et al.l (120081) . we find a mean difference 
(This Study - Literature) of -4.08 mA (a = 4.71 mA) (see 

4 T he results presented here are very similar to those used in IFriel et aLI 
1 2010!). The average difference in [X/Fe] is only 0.03 dex ± 0.01 dex (cr = 
0.06 dex). The largest differences were for Na (our [Na/Fe] values are lower 
by 0.14 dex) and Ca (our [Ca/Fe] values are lower by 0.18 dex). All other 
elements are within 0. 1 dex. 
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TABLE 6 

Equivalent Widths for Program Stars 



Wavelength (A) 


Species 


LEP (eV) 


logs/ 


Bel8 1163 


Bel8 1383 


Be21 T50 


Be21 T51 


Be22 414 


Be22 643 


Be32 16 


6300.30 


8.0 


0.00 


-9.72 




38.7 




38.9 


48.4 


114.0 


24.9 


6363.78 


8.0 


0.02 


-10.19 




20.8 




20.3 


27.3 


83.2 


16.3 


5688.19 


11.0 


2.11 


-0.42 














120.7 


6154.23 


11.0 


2.10 


-1.53 


60.8 


75.0 


74.6 


73.2 


77.8 




44.3 


6160.75 


11.0 


2.10 


-1.23 


83.8 


97.6 


97.9 


96.0 


103.6 




68.1 


5711.09 


12.0 


4.35 


-1.83 


122.5 












105.5 



REFERENCES. — Note. Table 6 is published in ils entirely in the electronic edition of the Astronomical Journal. A portion is shown here for guidance regarding its form and eontenl. 



TABLE 7 

Chemical Abundances for Berkeley 1 8 



TABLE 9 

Chemical Abundances for Berkeley 22 



Species 


Abundance 


(7 


N 


Abundance 


a 


N 


Species 


Abundance 


a 


N 


Abundance 


(7 


N 




Bel8 


-1163 




Be 18 


■1383 






Be22-414 




Be22-643 




[O/Fe] 








0.11 


0.07 


2 


[O/Fe] 


0.12 


0.04 


2 


0.21 


0.04 


2 


[Na/Fe] 


0.10 


0.05 


2 


0.32 


0.10 


3 


[Na/Fe] 


0.42 


0.23 


3 








rMa/Fel 


0.25 


0.10 


4 


0.10 


0.14 


4 


[Ma/Fel 


0.16 


0.12 


3 


0.19 


0.16 


4 


[Al/Fe] 


0.15 


0.04 


2 


0.22 


0.00 


2 


[Al/Fe] 


0.28 


0.00 


2 


0.49 


0.00 


2 


[Si/Fe] 


0.15 


0.18 


7 


0.19 


0.04 


10 


[Si/Fe] 


0.19 


0.07 


11 


0.19 


0.17 


5 


[Ca/Fe] 


-0.01 


0.18 


11 


0.25 


0.09 


6 


[Ca/Fe] 


0.23 


0.12 


6 


0.27 


0.27 


4 


[Ti/Fe] 


0.10 


0.17 


41 


0.27 


0.09 


30 


[Ti/Fe] 


0.31 


0.11 


17 


0.38 


0.19 


6 


[Mn/Fe] 


-0.27 


0.12 


3 


0.08 


0.08 


3 


[Mn/Fe] 


0.00 


0.10 


3 


-0.45 


0.09 


3 


[Fe I/H] 


-0.47 


0.10 


46 


-0.41 


0.12 


46 


[Fe I/H] 


-0.40 


0.11 


39 


-0.49 


0.22 


23 


[Fe II/H] 


-0.44 


0.17 


11 


-0.41 


0.18 


8 


[Fe II/H] 


-0.44 


0.09 


10 


-0.50 


0.33 


5 


[Co/Fe] 


0.12 


0.06 


3 


0.25 


0.08 


3 


[Co/Fe] 


0.25 


0.10 


3 


0.12 


0.03 


3 


[Ni/Fe] 


-0.07 


0.15 


6 


0.04 


0.12 


6 


[Ni/Fe] 


0.06 


0.16 


6 


0.14 


0.26 


6 


[Rb/Fe] 


-0.31 




1 


-0.15 




1 


[Rb/Fe] 


-0.20 




1 


-0.35 




1 


[Zr/Fe] 


-0.30 


0.09 


3 


-0.11 


0.03 


3 


[Zr/Fe] 


-0.14 


0.08 


3 








[Ba/Fe] 


0.25 




1 


0.34 




1 


[Ba/Fe] 


0.64 




1 


0.57 




1 


[La/Fe] 


0.28 


0.07 


2 


0.39 


0.00 


2 


[La/Fe] 


0.39 


0.07 


2 


0.36 


0.00 


2 


[Eu/Fe] 


0.25 




1 


0.34 




1 


[Eu/Fe] 


0.33 




1 


0.19 




1 
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TABLE 10 










Chemical Abundances for Berkeley 21 






Chemical Abundances for Berkeley 32 




Species 


Abundance 




N 


Abundance 


u 


N 


Species 


Abundance 


17 


N 


Abundance 


a 


N 




Be21-T50 




Be21-T51 






Be32 


■16 




Be32-18 




[O/Fe] 








0.18 


0.04 


2 


[O/Fe] 


0.15 


0.04 


2 


0.30 


0.07 


2 


[Na/Fe] 


0.39 


0.30 


3 


0.32 


0.18 


3 


[Na/Fe] 


0.22 


0.21 


3 


0.25 


0.28 


3 


[Mg/Fe] 


0.17 


0.13 


4 


0.20 


0.15 


4 


[Mg/Fe] 


0.22 


0.08 


4 


0.20 


0.08 


3 


[Al/Fe] 


0.15 


0.00 


2 


0.19 


0.07 


2 


[Al/Fe] 


0.16 


0.00 


2 


0.19 


0.04 


2 


[Si/Fe] 


0.18 


0.08 


10 


0.24 


0.08 


12 


[Si/Fe] 


0.18 


0.10 


13 


0.22 


0.05 


13 


[Ca/Fe] 


0.21 


0.12 


4 


0.19 


0.11 


6 


[Ca/Fe] 


0.16 


0.20 


1 1 


0.14 


0.12 


12 


[Ti/Fe] 


0.22 


0.05 


27 


0.17 


0.04 


30 


[Ti/Fe] 


0.23 


0.07 


26 


0.29 


0.08 


30 


[Mn/Fe] 


0.18 


0.08 


3 


0.03 


0.13 


3 


[Mn/Fe] 


-0.13 


0.15 


3 


0.03 


0.12 


3 


[Fe I/H] 


-0.26 


0.13 


44 


-0.35 


0.14 


31 


[Fe I/H] 


-0.38 


0.11 


59 


-0.37 


0.10 


52 


[Fe II/H] 


-0.26 


0.16 


12 


-0.33 


0.18 


6 


[Fe II/H] 


-0.42 


0.17 


14 


-0.36 


0.13 


16 


[Co/Fe] 


0.14 


0.08 


3 


0.15 


0.10 


3 


[Co/Fe] 


0.12 


0.10 


3 


0.18 


0.15 


3 


[Ni/Fe] 


0.07 


0.15 


7 


0.06 


0.15 


6 


[Ni/Fe] 


0.04 


0.12 


7 


0.03 


0.12 


7 


[Rb/Fe] 


0.05 




1 


-0.20 




1 


[Rb/Fe] 


-0.25 




1 


-0.10 




1 


[Zr/Fe] 


-0.05 


0.05 


3 


-0.16 


0.06 


3 


[Zr/Fe] 


-0.11 


0.03 


3 


-0.01 


0.03 


3 


[Ba/Fe] 


0.59 




1 


0.58 




1 


[Ba/Fe] 


0.37 




1 


0.20 




1 


[La/Fe] 


0.55 


0.07 


2 


0.57 


0.07 


2 


[La/Fe] 


0.42 


0.07 


2 


0.45 


0.07 


2 


[Eu/Fe] 


0.25 




1 


0.36 




1 


[Eu/Fe] 


0.23 




1 


0.38 




1 



Figure While the agreement is reason able, there is a 
small systematic difference. (We note that the Bragag lia et alj 
( 120081) spectra are of comparable spectral resolution, R = 
45,000, and similar S/N.) For Be32-18, uniformly increas- 
ing the EWs by 4.0 mA, while keeping all stellar param- 
eters unchanged, would increase [Fe/H] by 0.09 dex. We 
also compare our log g/ values for individual elements with 
iBragaglia et al.l ( 12008b (see Figure [T3ll. The largest difference 
is for Si, A logg/ = 0.12 dex. For all other elements, the av- 



erage differences in logg/ values for a given element are 0.06 
dex or smaller. 

Next, we compare our abundances [X/H] and [X/Fe] for a 
given cluster with literature values. In Figure [14] we compare 
our abundances for stars Be 21 T50 and Be 21 T51 w ith the 
three Be 21 stars analyzed by iHil l & Pasquini (1999). Our 
mean metallicity is [Fe/H] = -0.31 while their value is [Fe/H] 
= -0.54. For the elements in common, our [X/Fe] ratios are 
on average higher by 0.05 dex (a = 0.20 dex). For Ca, there is 



TABLE 11 
Chemical Abundances for 
PWM4 



Species 


Abundance 


(7 


N 


PWM4-RGB 1 


[O/Fe] 


0.18 




2 


[Na/Fe] 


0.42 


0.12 


3 


[Mg/Fe] 


0.23 


0.11 


4 


[Al/Fe] 


0.26 


0.04 


2 


[Si/Fe] 


0.10 


0.05 


7 


[Ca/Fe] 


0.19 


0.09 


4 


[Ti/Fe] 


0.14 


0.05 


6 


[Mn/Fe] 


-0.28 


0.13 


3 


[Fe I/H] 


-0.34 


0.15 


33 


[Fe II/H] 


-0.29 


0.12 


S 


[Co/Fe] 


0.13 


0.06 


3 


[Ni/Fe] 


-0.03 


0.11 


6 


[Rb/Fe] 


-0.30 




1 


[Zr/Fe] 


0.11 


0.06 


3 


[Ba/Fe] 


0.46 




1 


[La/Fe] 


0.22 


0.07 


2 


[Eu/Fe] 


0.12 




1 



TABLE 12 

Abundance Dependences on Model Parameters for Be 
21 T51 



Species 


r eff + 100K 


logg + 0.3 dex 


£, + 0.2 km s -1 


Total 


[O/H] 


-0.01 


0.11 


-0.01 


0.11 


[Na/H] 


0.09 


0.00 


-0.05 


0.10 


[Mg/H] 


0.03 


0.00 


-0.02 


0.04 


[Al/H] 


0.08 


-0.01 


-0.04 


0.09 


[Si/H] 


-0.04 


0.06 


-0.04 


0.08 


[Ca/H] 


0.10 


-0.02 


-0.10 


0.14 


[Ti/H] 


0.17 


0.01 


-0.10 


0.20 


[Mn/H] 


0.10 


0.02 


-0.17 


0.20 


[Fe I/H] 


0.06 


0.03 


-0.09 


0.11 


[Fe II/H] 


-0.09 


0.15 


-0.07 


0.19 


[Co/H] 


0.05 


0.06 


-0.06 


0.10 


[Ni/H] 


0.03 


0.08 


-0.07 


0.11 


[Rb/H] 


0.12 


0.00 


-0.02 


0.12 


[Zr/H] 


0.21 


0.01 


-0.03 


0.21 


[Ba/H] 


0.03 


0.07 


-0.17 


0.19 


[La/H] 


0.02 


0.12 


-0.06 


0.14 


[Eu/H] 


-0.01 


0.12 


-0.04 


0.13 


[O/Fe] 


-0.07 


0.08 


0.08 


0.13 


[Na/Fe] 


0.03 


-0.03 


0.04 


0.06 


[Mg/Fe] 


-0.03 


-0.03 


0.07 


0.08 


[Al/Fe] 


0.02 


-0.04 


0.05 


0.07 


[Si/Fe] 


-0.10 


0.03 


0.05 


0.12 


[Ca/Fe] 


0.04 


-0.04 


-0.01 


0.06 


[Ti/Fe] 


0.11 


-0.01 


-0.01 


0.11 


[Mn/Fe] 


0.04 


-0.01 


-0.08 


0.09 


[Co/Fe] 


-0.01 


0.03 


0.03 


0.04 


[Ni/Fe] 


-0.03 


0.05 


0.02 


0.06 


[Rb/Fe] 


0.06 


-0.03 


0.07 


0.10 


[Zr/Fe] 


0.15 


-0.01 


0.06 


0.16 


[Ba/Fe] 


-0.03 


0.05 


-0.08 


0.10 


[La/Fe] 


-0.04 


0.09 


0.03 


0.10 


[Eu/Fe] 


-0.07 


0.09 


0.05 


0.12 



a large difference between the two studies, A[Ca/Fe] = 0.30 
dex. For all other elements, the maximum difference in [X/Fe] 
is 0.11 dex. 

In Figure [15] we compare our Be 22 abu ndances for two 
stars w ith the two different stars analyzed by IVillanova et alj 
( 12005b . Our mean metallicity is [Fe/H] = -0.45 and their 
value is [Fe/H] = -0.32. For the elements in common, our 
[X/Fe] ratios are on average higher by 0.21 dex (a = 0.12 
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FIG. 12. — Comparis on of equivalent width measurements for Be 32 18 
between this study and Bragaglia et al. 1 2008). The one-to-one relation is 
shown (solid line) along with the linear fit to the data (dashed line). 
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FIG. 14. — Comparison of the abundance ratios [X /H] (upper) and [X/Fe] 
(lower) in Be 21 between this study (rilled circles) and Hill & Pasquini (1999) 
(crosses). The individual stars are marked and their relative positions are the 
same in each panel. 
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FIG. 15. — Comparison of the abundance ratios [X/ H] (upper) and [X/Fe] 
(lower) in Be 22 between this study (filled circles) and|Villanova et all <2005l) 
(crosses). The individual stars are marked and their relative positions are the 
same in each panel. 

In Figures [16] and [17] we compare our abundances for our 
two star s in Be32 with t h e mea n clu ster values from th e anal- 
yses of iBragaglia etal.1 (12008b and iFriel et all ( 120101) . Our 
mean metallici t y, [Fe /H] = -0.38, is 0.09 dex lower than 
IBragadia etail 120081) . [Fe/H] = -0.29. For elements in com- 
mon, our [X/Fel ratio s are on average higher than those of 
IBragadia etail ( 120081) by 0.06 dex (a = 0.07 dex). Our mean 



metall icity, [Fe/H] = -0.38, is 0.08 dex lower than IFriel et alj 
d2OT0b . [Fe/H] = -0.30. For elements in co mmon, our [X/Fe] 
ratios are on average higher than those of IFriel et al.l (120101) 
by 0.13 dex (a = 0.15 dex). This difference is driven primar- 
ily by Ti where A[Ti/Fe] = 0.43 dex. S uch a large abund ance 
difference is of concern, but as noted in lFriel et al.l (1201 Oh and 
apparent later in Figure [2T] [Ti/Fe] abundances appear to be 
especially vulnerable to at times large systematic differences 
between studies. For examp le, as we shall no te in Section 
4.4.1, M67 is in common to IFriel et ail d2010l) and Paper I, 
and while the two studies obtain very similar metallicities, 
A[Fe/H] = -0.01 dex, a large abundance difference exists for 
Ti, A [Ti/Fe] = 0.26 dex. The origin of these differences is not 
altogether clear, and does not appear to be due to simple ef- 
fects of a single source, but is likely driven by a complex com- 
bination of factors, such as EWs, stellar parameters, log gf 
values, specific lines utilized, reference abundances, analysis 
techniques, etc. 

In FigureQ~8] we comp are our abundan c es for the star Be32 
18 with the analysis of [Bragaglia et al.l ( 120081) . We derive 
a metallicity of [Fe/H] = —0.37 compared to their value of 
[Fe/H] = -0.27. The elements in common are Na, Mg, Al, Si, 
Ca, Ti, Ni, Ba, and Fe. For [X/H] we find a mean difference 
(This Study - Literature) of -0.05 dex ± 0.04 dex (a = 0.12 
dex). For [X/Fe] we find a mean difference (This Study - Lit- 
erature) of +0.05 dex ± 0. 04 dex (a = 0.13 d ex). We note that, 
as with the comparison to Fri el et ail (120101) . the agreement is 
especially poor for Ti, with A[Ti/Fe] = 0.27. 

From these comparisons, it is clear that there are system- 
atic differences at the ~0. 1 dex level between various stud- 
ies. For some elements, the abundance differences can be as 
large as 0.4 dex. Such differences should be borne in mind in 
the comparisons that follow. Of particular interest is that our 
metallicities are, with the exception of Be 21, lower than the 
literature values by -~0.1 dex which may potentially lead to 
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FIG. 17. — Same as Figure [T6l but for Ti, Fe, Co, Ni, Zr, and Ba. 
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FIG. 18. — Comparison of the abundance ratios [X/H] (up per) and [X/Fe] 
(lower ) in Be 32 18 between this study (filled circles) and IBragaglia et all 
<2008l) (crosses). 

larger [X/Fe] ratios than compared to the literature values. In 
the subsequent discussion, we will (at times) re-evaluate our 
conclusions by arbitrarily shifting our metallicities, [Fe/H], 
by +0.1 dex and decreasing our abundance ratios, [X/Fe], by 
0.1 dex. 

4. DISCUSSION 

4. 1 . Trends With Galactocentric Distance 

In Figure [19] we plot the metallicity, [Fe/H], versus Galac- 
tocentric distance. This figure includes measurements of 



open clusters from this paper (filled circles), Paper I (filled 
squares), and literature sources (p lus signs). (Altho ugh mem- 
bership of Be 3 1 remains unclear dFriel et al.ll2.Q10t) . we retain 
this key cluster, Rqc =12 kpc, in subsequent figures and dis- 
cussion.) The compilation of literature data (LTE abundances) 
and sources used to generate this figure can be found in Table 
[T3l The literature search was not exhaustive. We have not 
made a critical assessment of the literature sample and note 
that some of the abundance measurements were based on data 
that we would regard to have unacceptably low S/N. (In the 
Appendix A.l, we explore how selections based on S/N may 
affect the results.) This combined sample consists of some 
68 chemical a bundance measureme nts in 49 unique clusters. 
(We note that lPancino et al.l (120101) assembled a compilation 
of some 57 open clusters. They have more clusters in the so- 
lar vicinity but fewer in the outer disk, which is the focus of 
this series of papers.) One issue is that different authors may 
adopt independent distances. Theref ore, this figure shows 
the distances from ISalaris et al.l (120041) (upper panel) as well 
as the distances from the individual literature sources (lower 
panelfl In both panels, we join individual clusters with mul- 
tiple measurements using red lines. That is, in the event that 
a given cluster has been studied by two or more authors, we 
treat all studies as independent measurements regardless of 
the sampl e size, spectra l resolu tion, and S/N of each study. 
Following iPancino et al.l d2010l) . rather than averaging multi- 
ple measurements into a single value for a given cluster, we 
show all measurements to provide a "realistic idea of the un- 
certainties involved in the compilation." We also measure the 
lineaiQ fit to the local (Rgc < 13 kpc) and distant (Rcc > 
13 kpc) samples. The 13 kpc boundary dividing the local and 
distant clusters was arbitrarily chosen (see Appendix A. 2 for a 
more thorough examination of the division be tween the inner 
and outer disk and we als o refer the reader to I Jacobsonll2009l 
and lJacobson et al.l201 lb! for further discussion on this issue). 
As noted previously, the metallicity gradients in the two re- 
gions are considerably different. The inner region shows a 
steeper gradient than in the outer disk. Interestingly, the for- 
mal slopes, uncertainties, and dispersion about the linear fit 
are virtually identical in both panels, th at is, these va l ues do 
not depend upon the distance estimate, ISalaris et al.l (120041) 
vs. literature. Our conclusions would not change had we se- 
lected a boundary value of 1 2 kpc or 1 4 kpc rather than 1 3 kpc . 
For completeness, we note that the linear fit includes multiple 
measurements of a given cluster such that those clusters in 
effect are given more weight than a cluster with a single mea- 
surement. In subsequent figures, we use the literature distance 
estimates. 

The five clusters studied in the present pap er he on, or 
slightl y below, the mean trends. (If we adopt the lSalaris et al.l 
( 120041) distances, then the two most distant clusters in this 
paper are more metal-rich than the other three clusters. Al- 

5 When distances from Salaris et al. 12004) are not available, the values 
from the individual literature sources are used in the upper panel. Simi- 
larly, when distance s are not given in the literature sources, th e values from 
ISalaris et alj 120041) are employed. The ISalaris et alj (2004j) distances are 
based on main sequence fitting and are therefore sensitive to extinction and 
metallicity. The literature distances come from a variety of methods. As 
noted in Section 2.4, we use the red clump luminosity. Distances based on 
this approach are insensitive to metallicity, but identifying the red clump lo- 
cation and extinction affect our results. 

6 While we determine linear fits to the data throughout the paper, we are 
not suggesting that this is the correct function to use. Rather, we consider 
this a first pass to begin to understand the behavior of abundance ratios with 
distance, age, and/or metallicity, [Fe/H]. 
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TABLE 13 

Compilation of (LTE) Abundance Ratios for Open Clusters 



Age a Distance" Age b Distance 15 
Cluster Source (Gyr) (kpc) (Gyr) (kpc) [Fe/H] [O/Fe] [Na/Fe] [Mg/Fe] [Al/Fe] [Si/Fe] 



Be 17 


FriellO 


10.10 


11.70 


10.06 


10.89 


-0.10 


0.00 


0.33 


0.12 


0.25 


0.30 


Be 18 


This Study 


5.69 


13.10 


5.69 


12.09 


-0.44 


0.11 


0.21 


0.18 


0.19 


0.17 


Be 20 


Sestito08 


6.00 


16.37 


4.05 


16.12 


-0.30 




0.17 


0.28 


0.17 


0.05 


Be 20 


Y05 


4.10 


16.00 


4.05 


16.12 


-0.49 


0.18 


0.21 


0.24 


0.18 


0.06 


Be 21 


This Study 


2.18 


14.27 


2.18 


14.27 


-0.31 


0.18 


0.36 


0.19 


0.17 


0.21 



REFERENCES. — Bra gagliaOl = Bragaglia e t~al] 120011) : Bragaglia08 =[Bragaglia et al. t2008ft : Carraro04 = lCarraro etaD 120041) : Carrerall = 
ICarrera & Pancino i 2 0Tll); 

Carretta04 = |Carretta et al] {20041): Carretta05 =ICarretta et alj (20011); Carretta07 =|Carretta et aljf2007|): de Silva 06 = |De Silva etlri1f27j06l) : 

de Silva07 = |De Silva ct all <2007D: FriellO = |Friel et alj 120051) : I Jacobson et alj (2008, 2009): F riel et alj f201Cl) ; Hill99 = |Hill & Pasquinil fl999l) ; 

Jacobson 1 1 a = ( Jacobson et ajj |20JJj|) ; 

Jocobsonllb j Jacobson et aijj |2011b|) : Panc inolO =|Pan cino et al. i 20]]]); Paulson0 3 - Paulson e7al]{2003); Sestito08 = |Sestito et all f200l [20071 
[207)M Bi-agadia et al. i200j); Villanova05 =|VManova et al. (2005); Y05 = IYong et aljj2003) 



Note. Table 13 is published in its entirety in the electronic edition of the Astronomical Journal. A portion is shown here for guidance re- 
garding its form and content. 



a Ages and distances taken from the individual papers 

k Ages and distances taken fromfSalaris et al. i 2 0041) 

c Ages and distances taken from |Bragag lia&Tosilt200fl) 
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FIG. 19. — [Fe/H] vs. Galactocentric dist ance. In the upper panel, the 
distances are taken from Salaris et al. 1 2004]) while in the lower panel, the 
distances are from the literature papers. The data from this study are filled 
circles, the Paper I data are filled squares, and the literature data are plus 
signs. The connected red symbols represent multiple measurements of the 
same cluster. In both panels, we show the linear fit to the data with Rqq < 1 3 
kpc and Rqq > 13 kpc. A representative error bar is shown. 



though we note that this is not the case for our distance es- 
timates.) Given the inhomogeneous nature of this figure, it 
is difficult to ascertain whether or not the dispersion about 
the linear trend represents a real spread in metallicity or sys- 
tematic offsets between the various studies. Our typical mea- 
surement uncertainty for [Fe/H] is 0.15 dex, a value compa- 
rable with the dispersion about the linear fit which may in- 
dicate that the dispersion about the mean trend may be al- 
most entirely attributable to measurement uncertainties. In- 
deed, this issue of inhomogeneous comparisons is a recur- 
ring theme throughout the rest of the paper. We note that 



in this figure, and subsequent figures, no attempt has been 
made to normalize the abundances onto a common scale. It 
would be possible to make adjustments for the adopted so- 
lar abundances. However, quantifying systematic differences 
arising from equivalent widths, log gf values, stellar parame- 
ters, model atmospheres, spectrum synthesis codes, and so on 
would require an independent re-analysis of the entire data set 
(a substantial, but feasible analysis that should be conducted). 
Such a re-analysis, which shoul d also include non-LTE and/or 
3D effects when possible (e.g., iBergemann & Gehrenl 120081 : 
iLind et aT]|2009t) . is beyond the scope of the present paper. 

Had we shifted our metallicities (i.e., only our measure- 
ments in the present paper and in Paper I), [Fe/H], by +0.1 
dex, we would have obtained identical slopes for all subsam- 
ples considered in Figure [19] For the distant samples, the 
dispersion about the linear fit would decrease by 0.01 to 0.02 
dex depending on the adopted distance scale. 

Although we explicitly consider abundance trends with age 
in the following subsection, we now consider whether the 
trends and dispersions between metallicity and Galactocen- 
tric distance differ when considering open clusters of differ- 
ent age ranges. In Figuref20] we again plot metallicity versus 
Galactocentric distance, but the sample is split into various 
age groups: < 2 Gyr, 2 < Age < 5, and > 5 Gyr. These bins 
were arbitrarily chosen to contain roughly equal numbers of 
clusters and none of our conclusions would change had we 
changed the age groups by ± 0.5 Gyr. The upper panel of this 
figure shows the full sample along with the linear fit to the 
data for the local {Rgc < 13 kpc) and the distant (Rgc > 13 
kpc) samples. In the lower panels, each of which covers a dif- 
ferent age range, we superimpose the linear fit to the complete 
sample and we re-measure the dispersions about the linear fit 
for each subsample. The most notable feature of this dissec- 
tion is that for the local sample, Rgc < 13 kpc, the scatter 
about the mean relation is largest for the oldest clusters (> 5 
Gyr) and smallest for the youngest clusters (< 2 Gyr). This 
result would remain unchanged had we (i) shifted our metal- 
licities by +0. 1 dex or (ii) excluded the very old (10.2 Gyr) and 
metal -rich ([Fe/H] = +0.47) cluster NGC 6791 dCarretta et alj 
120071) . 

Clearly, it would be interesting to study whether the tran- 
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FIG. 20. — [Fe/H] vs. Galactocentric distance using the same symbols as in 
Figure [19] The upper panel shows all stars. The lower panels show stars in 
various age groups: < 2 Gyr, 2 < Age < 5, and > 5 Gyr. The linear fit to the 
data with Rqc < 13 kpc and Rqc > 13 kpc is shown in the upper panel. We 
superimpose these fits to the lower three panels but remeasure the dispersion 
about the linear fit. A representative error bar is shown. 

sition radius between the i nner and outer disk cha nges as a 
function of age. To this end. lJacobson et all (1201 lbl) explored 
this issue and found a suggestion of some variation. 

In Figure [21] we plot [X/Fe] for the a elements versus 
Galactocentric distance. In the bottom panel, a represents 
the direct average of O, Mg, Si, Ca, and Ti. As in previous 
figures, we compute and plot the fit to the data for Rgc < 13 
kpc and Rgc > 13 kpc. For all elements, there is no statisti- 
cally significant difference between the gradient for the local 
and distant samples. While combining the a elements into a 
single measure may be a convenient way to examine their be- 
havior and to facilitate comparisons with model predictions, 
we note that the individual elements do not appear to follow 
identical patterns: [Ca/Fe] is roughly solar at all distances, 
[Si/Fe] and [Mg/Fe] are uniformly enhanced at all distances, 
and [O/Fe] and [Ti/Fe] show considerably higher ratios in the 
outer disk compared to the inner disk. Thus the behavior of 
[a/Fe] would appear to be driven primarily by O and Ti. As 
noted in Paper I, the outer disk open clusters are probably not 
members of an old er population such as the thick disk (e.g., 
iBensbv e"ta l. 2004) despite some similarities in [X/Fe] abun- 
dance ratios. 

We plot the light elements (Na and Al; Figure l22l . the 
Fe-peak elements (Mn, Co, and Ni; Figure l23l . and the 
neutron-capture elements (Zr, Ba, La, and Eu; Figure l24l 
versus distance. For the neutron-capture elements discussed 
here, we note that Zr, Ba, and La are predominantly pro- 
duced via the ^-proces s while Eu is produced via the r-process 
(ISimmerer et al. 1 120041) . For the neutron-captur e elements Zr, 
Ba, an d La, we include the abunda nce data from lD'Orazi et all 
( 120091) and iMaiorca et al.l (120111) and dis t ances, when not 
listed in those re ferences or ISalaris et al.l (120041) . from the 
IChen et all (120031) catalogue. We again compute and over- 
plot the linear fit to the local (Rgc < 13 kpc) and distant (Rgc 
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FIG. 22. — Same as Figure |2~T1 but for the light elements [Na/Fe] and 
[Al/Fe] . 



> 13 kpc) samples. While most elements show no significant 
difference in gradient between the local and distant samples, 
within the limited data there is a suggestion that Mn, Ba, La, 
and Eu have different gradients between the local and distant 
samples. 

For all a elements, the gradients are positive (or zero). The 
maximum value, 0.02 ± 0.01 dex kpc -1 , occurs for O and Si 
in the range Rgc < 13 kpc. Therefore, the a elements show 
only a small trend with increasing Rgc- Similarly, Na and 
Al have positive (or zero) gradients suggesting a slight trend 
with Rgc- The maximum value is 0.03 ± 0.01 dex kpc" 1 for 
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FIG. 23. — Same as Figure |2T1 but for the Fe-peak elements [Mn/Fe], 
[Co/Fe], and[Ni/Fe]. 
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FIG. 24. — Same as Figure l2"T"l but for the neutron-capture elements [Zr/Fe], 
[Ba/Fe], [La/Fe], and [Eu/Fe]. 

Al in the range Rgc < 13 kpc. For the Fe-peak elements Mn, 
Co, and Ni, the gradients are both positive and negative. The 
largest value is for Mn in the range Rgc < 13 kpc, -0.06 ± 
0.01 dex kpc" 1 . All other gradients are small. The neutron- 
capture elements Ba, La, and Eu show large positive gradients 
in the range Rqc < 13 kpc. In the outer disk, Rgc > 13 kpc, 
their gradients are relatively flat with maximum values of 0.02 
dex kpc -1 . Zr shows no trend with Galactocentric distance. 

Another way to compare the abundances between the inner 
and outer disk is to consider the mean [X/Fe] ratios in the two 
regions. We again choose a boundary value of Rgc = 13 kpc 
between the inner and outer disks. We then measure the mean 
[X/Fe] ratio as well as the standard error of the mean. By 
adding the two standard errors of the mean in quadrature and 
comparing this value to the difference between the two means, 
we obtain the level of significance. In Table [T4l we show for 
each element the level of significance that the [X/Fe] ratios 
differ between the inner and outer disk. The elements that 



differ at the 3-rj level or higher are O, Na, Al, and Ti. While 
we note that a is also significant at the 3-er level, Mg and Si 
are less than 1-er and Ca is less than 2-a. In this table, we 
also show results when using boundary values of 10 kpc and 
15 kpc. For some elements, the level of significance depends 
strongly upon the adopted boundary. 

We now re-consider our results when we arbitrarily in- 
crease our [Fe/H] values by 0.1 dex, and therefore decrease 
our [X/Fe] by 0. 1 dex. Recall that this shift in [Fe/H] is moti- 
vated by differences between our [Fe/H] values and literature 
values for clusters in common with the literature. We find 
that while there remain differences in the mean abundances 
between the inner and outer disks, the level of significance is 
lower in some cases; O, Na, Al, Zr, and a are significant at 
the 2-a level and Ti is significant at the 4-er level. 

4.2. Trends With Age 

Next we explore whether there are any trends between 
metallicity, [Fe/H], and abundance ratios, [X/Fe], versus age 
(Figures [25] and |26"1 ). In these figures we adopt the ages 
from lSalaris et al] " il2004l) and individual clusters with multi- 
ple measurements are connected wi t h red lines. For the open 
clusters not included in lSalaris et al.l (120041) . we adopt the ages 
given in the literature sources (see Table [T3l for all ages and 
references). In all panels of these figures, we compute the 
linear fit to the data and show the gradient, uncertainty, and 
dispersion about the linear fit. For all elements there are no 
significant trends, at the >3-a level, between abundance and 
age. The most significant trends we found were for Co, +0.03 
± 0.01 dex/Gyr and Ba, -0.03 ± 0.01 dex/Gyr. (Lines from 
both elements are affected by hyperfine and/or isotopic split- 
ting and therefore these elements are more likely to be af- 
fe cted by systematic un certainties than other elements.) 

iD'Orazi et"aj] (120091) found that the [Ba /Fe] ratio decreased 
with increasing age. Maj orca et al.l (1201 ll) extended the mea- 
surements to additional neutron-capture elements and con- 
firmed the trends between decreasing [X/Fe] with increasing 
age. We find a similar trend, albeit with a smaller amplitude 
and significance for Zr and Ba. A likely explanation for this 
result is that our new measurements of neutron-capture ele- 
ments in our sample of older clusters fall a t hi gher [X/Fel ra- 
tios th an the values in lD'Orazi et alJ (120091) and lMaiorca et al] 
d201 lb . For La, we find an opposite trend of increasing 
[La/Fe] with increasing age. 

For some elements, primarily Ni, the dispersion about the 
linear fit {a = 0.07 dex) is comparable to the measurement 
uncertainty (0.06 dex). We also measure the mean abundance 
ratio and dispersion for the three age groups used earlier: < 
2 Gyr, 2 Gyr < Age < 5 Gyr, and > 5 Gyr. Given the lim- 
ited data and inhomogeneous analyses, it would be prema- 
ture to conclude that the dispersion differs between the vari- 
ous age subsamples, but there are hints for some elements of 
unusually large, or small, dispersions in various age subsam- 
ples. The presence, or absence, of such trends and dispersions 
would be revealed from a homogeneous analysis. 

Again, we reconsider the results when we arbitrarily in- 
crease our metallicities by 0. 1 dex, and decreasing our [X/Fe] 
ratios by 0. 1 dex. For all elements, the gradients change by 
<0.01 dex/Gyr and the dispersion about the linear trends gen- 
erally decreases by a very small amount. 

4.3. Trends With Metallicity 

In Figure [27] we plot the abundance ratios [X/Fe] versus 
metallicity, [Fe/H]. In all panels, we plot the linear fit to the 
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TABLE 14 

Mean [X/Fe] Ratios, Standard Error of the Mean, and Level of Significance When Comparing the Inner and Outer Disk For 

Different Boundary Values. 
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FIG. 25. — Abundan ce r atios [Fe/H] and [X/Fe] vs. age. The symbols are 
the same as in Figure \\9\ The dashed blue line shows the linear fit to the 
data. In each panel, we also show the mean abundance ratio and dispersion 
for three age groups: (1) < 2 Gyr, (2) 2 < Age < 5, and (3) > 5 Gyr. A 
representative error bar is included in each panel. 



data and show the slope, uncertainty, and dispersion about the 
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are given more weight than a cluster with a single measure- 
ment.) W e overplot solar ne i ghborh ood thin and thick disk gi- 
ants from lAlves-Brito et alj ( 120101) as well as solar neighbor- 
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FIG. 26. — Same as Figuref25]but for different elements. 
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hood giants from iLuck & Heiterl d2007) which are primarily 
thin disk objects. The open cluster abundance measurements 
come almost exclusively from giant stars. Therefore, when 
comparing the clusters with field stars, we chose to include 
only giants in order to minimize any systematic abundance 
differences. For example, comparison of the results for Mg 
betwe en bu lge and local thick disk s tars from lFulbright et all 
(2007|) and|Ai ves-Brito et al.l d2010l) illustrates the potential 
systematic differences that may arise when comparing dwarfs 
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with giants. Additionally, elements such as Na have impor- 
tant non- LTE effects whic h differ between dwarfs and gi- 
ants (e.g., iLind et all 1201 lb . Arguably, any dwarf versus gi- 
ant abundance discrepancies may be small relative to the sys- 
tematic differences arising from an inhomogeneous compar- 
ison (stellar parameters, atomic data, equivalent widths, so- 
lar abundances, model atmospheres, spectrum synthesis soft- 
ware, etc). 

For almost all elements, the trends between [X/Fe] ver- 
sus [Fe/H] for open clusters and field stars are in fair to 
good agreement. Such agreement would suggest, at face 
value, that the open clusters, which span a range in ages 
and distances, have experienced a similar chemical evolution 
history as did the local field stars which span an unknown 
(but presumably large) range in age. The obvious discrep- 
ancies include Mn and Ba (these elements are affected by 
hyperfine an d/or isotopic sp l itting) and perhaps O and Na. 
For Mn, the iLuck & He"iteri (I20Q7T) field giants have higher 
[Mn/Fe] ratios than do the open clusters. For Ba, the field 
giants hav e lower [Ba/Fe ] ratio s than do the open clusters. 
As noted, iD'Orazi et al.l (120091) found a trend of decreas- 
ing [Ba/Fe] with increasing age in the open clusters and the 
lower [Ba/Fe] ratios in field stars could arise if the field stars 
are systematically older than the open clusters. In the open 
clusters, [O/Fe] appears lower than field stars of comparable 
metallicity while [Na/Fe] appears systematically higher than 
in field stars. While such abundance patterns resemble the 
O-Na anticorrelation seen in globular clusters dGratton et alJ 
l200l.the open clust ers do not display this abundance pattern 
(Ide Silvaet ai] l2009). The open cluster data are inhomoge- 
neous as are the comparison field stars. Until a systematic and 
homogeneous analysis is performed upon a sample of open 
clusters and field stars, we cannot definitively say whether or 
not the open clusters and field stars have the same, or differ- 
ing, chemical abundance trends. 

In Figure [28] we compare the predicted error in [X/Fe] due 
to uncertainties in the model parameters (i.e., the values in the 
fourth column of Table [T2l i with the dispersion about the lin- 
ear fit to open clusters in the [X/Fe] versus [Fe/H] plan<£] For 
O, Si, Ti, Ni, and Eu, these values are in good agreement sug- 
gesting that the dispersion about the linear fit to [X/Fe] versus 
[Fe/H] can be entirely attributable to the measurement uncer- 
tainties. Given the inhomogeneous nature of the open clus- 
ter data, this is a slightly surprising result. Notably, [Ni/Fe] 
and [a/Fe] show dispersions about the linear fit of only 0.07 
and 0.06 dex, respectively. For other elements, notably Co, 
Zr, Ba, and La, the predicted error in [X/Fe] is considerably 
smaller than the measured dispersion (which reaches values as 
high as 0.24 dex for Zr and Ba). Such a result would indicate 
that these elements show a real dispersion amongst the open 
clusters or that the combined sample is significantly affected 
by systematic offsets between the various individual studies. 
For Ba, the lines are generally quite strong and require con- 
sideration of hyperfine splitting such that the dispersion will 
likely be larger than the errors arising solely from uncertain- 
ties in the stellar paramete rs. Furthermore, the t rend between 
[Ba/Fe] and age found by ID'Orazi et al.l (120091) would intro- 
duce additional scatter when considering clusters of all ages 
in the [Ba/Fe] vs. [Fe/H] plane. Finally, it becomes more dif- 

7 This approach is only meaningful if we expect the dependence of [X/Fe] 
versus [Fe/H] to be linear. Inspection of Figure l27l would suggest that the dis- 
tributions are often non-linear. In addition, we need to be mindful of potential 
systematic differences between various studies which could serve to increase 
the observed element abundance dispersion. 



ficult to accurately measure the strengths of strong lines such 
as Ba and their strength also means that the lines form at shal- 
lower optical depths such that the LTE assumption is less re- 
liable. 

La may offer a more reliable measure of the s-process since 
there are more lines available to measure including a number 
of weaker lines. L a may show a trend with age (a negative 
trend according to dMaiorca et al.ll201l1) and a positive trend 
according to this study) and therefore inclusion of clusters 
with a large range of ages may serve to increase the dispersion 
for [La/Fe] vs. [Fe/H] beyond that expected only from uncer- 
tainties in the stellar parameters. Finally, the large dispersions 
for Zr and Co may be due, in part, to the small number of lines 
typically used to derive abundances (although some other ele- 
ments also measured from comparably small numbers of lines 
show agreement between the predicted and observed abun- 
dance dispersions). 

We re-examine our results when increasing [Fe/H] by 0. 1 
dex and decreasing [X/Fe] by 0.1 dex. For the linear fit to 
the [X/Fe] vs. [Fe/H] data, the slopes increase by an average 
of 0.08 ± 0.01 dex/dex (i.e., the negative slopes become flat- 
ter). Such a result is readily anticipated given that our clus- 
ters are amongst the most metal-poor in the combined sam- 
ple and lowering the [X/Fe] ratios will therefore produce flat- 
ter slopes when fitting [X/Fe] vs. [Fe/H]. Regarding Ti, our 
highest [Ti/Fe] ratios are comparable to the upper envelope of 
field star data and some open cluster literature [Ti/Fe] data lie 
substantially below the relation defined by the field star data. 
When adjusting our [Fe/H] and [X/Fe] ratios by 0.1 dex, the 
dispersion about the linear fit is essentially unchanged and the 
clusters and field stars remain in good agreement with the ex- 
ception of the elements already noted, Mn and Ba and perhaps 
O and Na. 

4.4. The Formation and Evolution of the Outer Galactic 

Disk 

AAA. Insights from Open Clusters 

We now comment on what we can learn about the origin and 
evolution of the outer disk based on metallicities and chemical 
abundance ratios in open clusters. The existing data indicate 
that the outer Galactic disk open clusters (Rgc > 13 kpc) are 
uniformly metal-poor, i.e., [Fe/H] < 0.0, albeit with metal- 
licities higher than expected based on an extrapolation of the 
abundance gradient of local open clusters (Rgc <13 kpc). In- 
terestingly, we note that the handful of fi eld red giants in the 
inner Galactic disk dBensbv et alj|2010l) are also metal-poor 
and do not lie on a linear extrapolation of the metallicity gra- 
dient based on the solar neighborhood, although these objects 
have unknown ages and may belong to the thick disk rather 
than the thin disk. All recent studies of open clusters confirm 
that the metallicity gradient in the outer disk is flatter than the 
metallicity gradient near the solar neighborhood. As noted 
in the introduction, some external galaxies exhibit similar be- 
havior in their outer disks. 

We also note that almost all outer disk open clusters have 
super-solar [a/Fe] and [Eu/Fe] ratios. Such abundance ra- 
tios are different from the values found in the more metal-rich 
open clusters in the solar neighborhood. As already discussed, 
the [O/Fe], [Ti/Fe], and [a/Fe] ratios differ between the inner 
and outer disk at the 3-er level or higher (assuming a boundary 
of 13 kpc). Mg, Si, and Ca do not show significant differences 
for [X/Fe] between the two regions. Depending on the choice 
of the boundary value, [Eu/Fe] differs between the inner and 
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FIG. 27. — [X/Fe] vs. [Fe/H]. Literature open clusters (plus signs), clusters from this study (cl osed black circles), and clusters from Paper I (closed black 
squa res) are plotted. We ov erplot thin disk (black stars) and thick disk (blue triangles) giants from Alves-Brito et al. (2010) as well as (mainly thin disk) giants 
fromlLuck & Heiter 1 2007) (red open circles). A representative error bar is shown. We plot the linear fit to the open cluster data and write the slope, uncertainty, 
and dispersion about the slope. 



outer regions at the 5.0-cr, 2.0-cr, and 0.9-cr level for 10 kpc, 
13 kpc, and 15 kpc, respectively. 

As noted, the possibility exists that our [X/Fe] measure- 
ments are systematically high. In order to produce [X/Fe] 
ratios for O, Ti, and Eu that agree within 2-<r between the in- 
ner and outer disk would require decreases of 0.2, 0.4, and 0. 1 
dex for [O/Fe], [Ti/Fe], and [Eu/Fe], respectively. Additional 
data and a homogeneous re-analysis of the existing data may 
be necessary to clarify this intriguing situation of enhanced 
[a/Fe] and [Eu/Fe] ratios in the outer disk. 

All 5-process elements studied show a significant abun- 
dance scatter in the outer disk (and in the solar neighbor- 
hood) extending to values as high as [X/Fe] ~ +1.0. Such 
high abundances of the s-process elements would suggest that 
asymptotic giant branch stars have played, in some cases, a 
major role in the chemi cal evolution of open clusters relative 
to field stars. Recently. iMaiorca et al.1 (120121) have suggested 
that the large [s-element/Fe] ratios seen in young metal-rich 



open clusters can be explained by low-mass AGB stars with 
efficient s-processing. 

When comparing the distant open clusters with solar neigh- 
borhood field stars at a given metallicity, the abundance ratios 
[a/Fe] and [Eu/Fe] are in fair agreement. Such similarity in 
abundance ratios would suggest that in the range of metallici- 
ties spanned by both sets of objects, the relative contributions 
from Type II supernovae (SNe II) and Type la supernovae 
(SNe la) are similar. The primary difference is that the chem- 
ical enrichment in the outer disk did not yet reach the metal- 
licities of the solar neighborhood. Therefore, the conclusion 
we drew in Paper I that "high abundances of a-elements in- 
dicate rapid star formation, such that Type la supernovae did 
not have sufficient time to evolve and contribute to the chem- 
ical evolution" may not be correct since it now appears likely 
that SNe la have made a contribution. Two differences be- 
tween the present work and that of Paper I are (1) the number 
of open clusters with chemical abundance measurements has 
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FIG. 28. — Predicted error in [X/Fe] due to uncertainties in the model pa- 
rameters, <7 (model), v ersus the dispersions about the linear fit to the open 
cluster data from Figure l27l a (gauss). The one-to-one relation is shown. 

increased considerably and (2) the comparison field stars are 
all giants which may remove any systematic abundance differ- 
ences arising from dwarf versus giant comparisons. Clearly it 
would be of great interest to identify and analyze more metal- 
rich distant open clusters. 

The metal-poor open clusters in the outer disk have super- 
solar ratios of the a elements and Eu which are primarily pro- 
duced in massive stars. Since we do not see significant con- 
tributions from SNe la, we may conclude that the star forma- 
tion in the outer disk was not prolonged otherwise we would 
have seen [a/Fe] substantially lower than in field stars of com- 
parable metallicity as found in the more metal-rich star s of 
the nearby dwarf spheroidal galaxies dTblstov et al.l2009l) . As 
discussed in Paper I, we would naively expect that the lower 
density in the outer disk would result in a slower star forma- 
tion rate relative to the solar neighborhood. In this scenario, 
we would expect the metallicity gradient to be continuous and 
that at a given metallicity, the outer disk would have lower 
[a/Fe] ratios than in the solar neighborhood. As we spec- 
ulated in Paper I, one possibility for these somewhat unex- 
pected characteristics in the outer disk is that a merger event 
and/or infall of material triggered a burst of star formation. 
Such a possibility was reinforced by the results from Paper 
III in which we found that the young outer disk Cepheids are 
more metal-poor than the older outer disk open clusters and 
that some of the Cepheids had enhanced [a/Fe] ratios. 

Flat abundanc e gradients could also be produced by ra- 
dial m i xing (e.g.. | Roskar et al.l 2008; Sa nchez-Blazquez et al.l 
I2009t iMinchev et al.l 1201 U 120121) and it is suspected 
that these proc esses may play a role in producing the 
thick disk (e.g.. ISellwood & Binnevll2002t iHay wood! l2008t 



ISchonrich & Binney||2009l) . Uflkova et al.l (120121) have shown 
that the Galactic bar and spiral arms may be responsible for 
moving NGC 6791 from the inner disk to its current location, 
although they regard the probability to be very low. While 
it may be possible for strongly bound clusters to survive the 
migration, to our knowledge no detailed study has yet been 
performed. 

Finally, we caution that any interpretation of the abundance 



ratios from compilations needs to acknowledge that the open 
clusters were studied by various authors who adopt different 
analysis techniques. We also note that for M67, a recent study 
bv lOnehag et al.l (1201 ll) showed that this cluster has a metal- 
licity [Fe/H] = +0.02 and abundance ratios [X/Fe] within 0.03 
dex of the solar values. Their strictly differential analysis of 
one dwarf star with almost identical stellar parameters to the 
Sun, a so-called solar twin, enabled very high precision. For 
comparison, in Paper I our analysis of M67 giants found near 
solar metallicity, [Fe/H] = 0.02, but [X/Fe] ratios that differed 
from the solar value. The elements with [X/Fe] > 0.1 dex 
are Na, Mg, Al, Ti, and La and those with [X/Fe] < -0.1 
dex are Mn and Zr. Table [T3l includes an additional study 
of M67 giants by Friel & Jacobson (e.g.. | Friel et al.l I2005t 
I Jacobson et all 120081 120091: iFriel et alJEOlOt Uacobson et al l 
2011a) that also found [X/Fe] ratios that differ from solar 
by 0. 1 dex or more, including Ti for which they find [Ti/Fe] 
< -0. 1 . Therefore, the lOnehag et al.l d201 ll) results serve to 
highlight any systematic abundance offsets between studies 
that have included this target. We also note that our Paper 
I analysis included the bright giant Arcturus such that future 
studies can identify abundance offsets. Additionally, mem- 
bership has not been unambiguously d etermined for some key 
clusters, e.g., Be 31 (IFriel et al.ll2010l) . 

4.4.2. Insights from Field Red Giant Stars 

Chemical abundance measurements exist for a number 
of field stars beyond Rqc = 10 kpc. These objects in- 
clude Cep heids (which we shall d iscuss in Sec 4.4.3), young 
OB stars fDaflon & Cunhal 120041) as well as field red giants 
dCarnev et al.l l2005bt iBensbv et alj|201 ll) which we assume 
to have ages comparable to the open clusters. We focus here 
on the field red giants and compare the abundances obtained 
in these objects to the open cluster data. The main issue we 
seek to address here is whether the field star population is 
chemically distinct from the open clusters. While we do not 
necessary expect any chemical differences between the two 
samples, it is important to conduct the comparison. 

Before continuing, we caution that outer disk field red gi- 
ants may be affected by selection biases. Stars beyond say 
^gc =10 kpc are likely chosen to lie well above, or below, 
the Galactic plane to avoid reddening. We speculate that such 
selection criteria may therefore give strong w eight towards 
thick disk stars over thin disk stars (although see lBensbv et al.l 
1201 ll for a discussion of the scale length and scale height of 
the thick disk). 

In Figure |29] we plot [Fe/H], [Mg/Fe], [Si/Fe], and 
[Ti/Fe] vs. Rqc (left panels) for open cl usters and field stars 
Carnev et all2005bl:lBensbv et alj|201 ll) . In these panels, we 
plot the linear fit to the local and distant open cluster samples. 
The first nota ble aspect from this f igur e is that the field star 
samples from lCarney et al.l (12005b!) and lBensby et all (120111) 
do not extend to the same large distances as do the open clus- 
ters. The second point we highlight is that the field stars 
are, on average, more metal-poor than the open clusters at 
the same distances. For open clusters with 10 kpc < Rqc < 
12 kpc the mean metallicity is [Fe/H] = -0.20 and the disper- 
sion is cr[Fe/H] = 0.11 dex. For the field stars in the same 
range of Galactocentric distances, we find a mean metallicity 
of [Fe/H] = -0.48 and dispersion of er[Fe/H] = 0.13 dex. It 
is not obvious why the metallicity difference exists between 
the open clusters and field stars at the same Galactocentric 
distances. Systematic offsets between the various studies at 
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FIG. 29.— [Fe/H], [Mg/Fe], [Si/Fe], and [Ti/Fe] vs. R G c (left panels) and 
[Mg/Fe], [Si/Fe], and [Ti/Fe] vs. [Fe/H] (rig ht panels) for open clusters (plus 
signs ), field stars fr om Carney et al. ( 2005b) (filled red circles), and field stars 
from Bensby et al. 1 201 1) (blue open triangles). In the left panels, we fit the 
local (Rac < 13 kpc) and distant (Rqc > 13 kpc) open clusters and in the right 
panels, we fit the open cluster data. Representative error bars are included in 
each panel. 



the ^0.3 dex level would be necessary to explain these dif- 
ferences or the field stars could be thick disk objects. Alter- 
natively, it may be that the difference in metallicity between 
the samples are due to selection biases in the field stars. Rais- 
ing our open cluster metallicities by 0. 1 dex would only serve 
to increase the discrepancy in [Fe/H] between open clusters 
and field stars. The third poin t we no te is that the dispersion 
in [Fe/H] for the lBensbv et al.l (1201 11) data (0.13 dex) is com- 
parable to the dispersion in [Fe/H] for the open clusters at 
similar distances (0.11 dex). 

We now compare the distributions of [Mg/Fe], [Si/Fe], and 
[Ti/Fe] vs. Rcc for the field stars and open clusters. Taking 
objects in the range 10 kpc < Rcc < 12 kpc, we find mean val- 
ues and dispersions of [Mg/Fe] = 0.15 (er = 0.09) for the open 
clusters and [Mg/Fe] = 0.11 (<r = 0.09) for the field stars. For 
Si, we find mean values and dispersions of [Si/Fe] = 0.16 {a 
= 0.08) and [Si/Fe] = 0.12 (a = 0.06) for the open clusters and 
field stars, respectively. For Ti, we find mean values and dis- 
persions of [Ti/Fe] = -0.03 (a = 0.15) and [Ti/Fe] = 0.09 (a = 
0.09) for the open clusters and field stars, respectively. While 
[X/Fe] is in very good agreement between open clusters and 
field stars in the same range of Galactocentric distances, there 
are very few field stars beyond Rcc =13 kpc with chemical 
abundance measurements. There is a clear need to study more 
distant field stars to understand the chemical properties of the 



outer Galactic disk. Had we decreased our [X/Fe] ratios by 
0. 1 dex, arguably the agreement in [X/Fe] between field red 
giants and open clusters would remain satisfactory. 

In the right panels of Figure [29] we compare [X/Fe] vs. 
[Fe/H] between the field stars and the open clusters. It seems 
that the field stars have lower [X/Fe] ratios, on average, than 
do the open clusters at the same metallicity. Within the limited 
samples and bearing in mind the measurement uncertainties 
and possible systematic offsets, the field stars and open clus- 
ters seem to follow the same trends between abundance ratios 
[Fe/H] and [X/Fe] vs. R GC and between [X/Fe] vs. [Fe/H]. 
Therefore, our tentative conclusion is that the interpretation 
of the outer Galactic disk does not depend significantly upon 
whether we use open clusters or field stars although we re- 
mind the reader that selection biases for the field stars may 
give much stronger weight to thick disk objects. 

4.4.3. Insights from Cepheids 

Studies of the chemical abundances of Cepheids provide 
a complementary view on the radial abundance gradients in 
the disk. With their high masses and short lifetimes, the 
Cepheids likely reflect the present-day chemical composition 
of the disk. A number of studies within the last decade have 
explored the radial abundance gradients as traced by Cepheids 
(e.g..lAndrievsky et al.l[2002b1faE l2004t jLemasle et alj|2008t 



Luck et al.l2003l.l201 ltlEuck & Lamberfeoi UlPedicelli et al.l 



2009tlYong et al.ll2006l) . 

While the time variation in radial abundance gradients can 
be obtained by comparing the younger and older open clus- 
ters, we note that the sample sizes are modest (e.g., our com- 
pilation has 24 and 43 clusters with ages below 1 Gyr and 
above 2.5 Gyr, respectively, and of the 24 clusters younger 
than 1 Gyr, only 3 lie beyond Rcc = 12 kpc and none beyond 
14 kpc). Instead, comparison between Cepheids and old open 
clusters offers a larger sample (e.g., iLuck & Lambert! [201 ll 
have a homogeneous sample of 339 Cepheids including 40 
and 15 objects beyond 12 kpc and 14 kpc, respectively) and 
arguably a more robust measure of the time variation of ra- 
dial abundance gradients. (We note that planetary nebulae 
also offer important insights into the chemical evolution of 
the Galactic disk, although their distant uncertainties are large 
relative to those of open clusters and Cepheids. We refer th e 
inte rested reader to iMaciel et al.l (120031) . ICosta et al.l (120041) . 
an d lMaciel& Costal (120091) and references therein.) 

One concern is that there may be large systematic differ- 
ences between the abundances obtained from open cluster 
red giants and those from Cepheids. Indeed, caution must 
be exercised when directly comparing abundances between 
open clusters and Cepheids. However, abundance gradients 
basically represent a differential abundance comparison be- 
tween similar stars but at different Galactocentric distances. 
If each sample (open cluster giants or Cepheids) is analyzed 
uniformly, then the systematic errors in the analysis of a par- 
ticular type of star should largely cancel thereby enabling the 
radial gradient comparisons we seek to conduct. 

In Figure [30] we plot [Fe/ H] vs. Rcc for the open clus- 
ters and for the Cepheids from lLuck & Lambert! d201 ll) . This 
figure enables us to address whether the older open clus- 
ters have a different [Fe/H] radial gradient than the young 
Cepheids. Therefore, we only plot open clusters with ages 
greater than 2.5 Gyr. Additionally, we divide the open clus- 
ters and Cepheids into a local (Rcc < 13 kpc) and distant 
(Rgc > 13 kpc) samples. We measure the linear fit to all 
subsamples and write the slope, uncertainty, and dispersion 
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FIG. 30. — [Fe/H] vs. Rqq for the older open clusters (age > 2.5 Gyr; left 
panel) and Cepheids (right panel). Data with Rqq < 13 kpc are in grey while 
data with Rqc > 13 kpc are in red. Open clusters from our work are show n 
as filled circles while the Cepheid data are from Luck & Lambert 12011). 
We determine the linear fit, uncertainty, and dispersion about the linear fit to 
the local and distant open clusters and Cepheids. For the Cepheid data, we 
superimpose the linear fits to the open cluster data. 



about the linear fit in Figure [30] As noted in previous stud- 
ies of Cepheids, there is no sudden change in gradient be- 
tween [Fe/H] and Rgc- In contrast to the open clusters, the 
Cepheids decrease smoothly in metallicity with increasing 
distance. The same result would apply had we overplotted 
our Cepheid da ta from Paper III. Since 21 of those 24 stars 
were studied bv lLuck & Lambertl (1201 lb . for convenience we 
rely only upon their data to ensure a homogeneous sample. 
Next, we note that for the inner samples, the older open clus- 
ters show a much steeper abundance gradient (-0.12 ± 0.01 
dex/kpc) than the younger Cepheids (-0.06 ± 0.00 dex/kpc). 

In Figure [31] we plot [Mg/H] and [Mg/Fe] vs. Galacto- 
centric distance as well as [Mg/Fe] vs. [Fe/H] for the open 
cluster and Cepheid samples. We again determine the lin- 
ear fit to the local and distant samples. For the local sample 
(Rgc < 13 kpc), the older open clusters show a steeper gra- 
dient for [Mg/H] compared to the younger Cepheids, a result 
also seen for [Fe/H]. This analysis was repeated for each el- 
ement and in Figure [32] we compare the gradients for [X/H] 
vs. Rqc between open clusters and Cepheids with Rgc < 13 
kpc. For the open clusters, we consider only clusters with 
ages > 2.5 Gyr (upper panel) and > 1.5 Gyr (middle panel). 
For both cases, we find that the gradient for [X/Fe] vs. Rgc 
is always shallower for the Cepheids compared to the older 
open clusters, with the exception of La. Such a time variation 
in radial abundance gradient is in th e same sense as predicted 
by chemical evolution mode ls (e.g.. lHou et al.ll2000t iFu et"ai1 
l2009t IPilkingto n et al.ll2012l) and is believed to be due to the 
"inside-out" formation of the disk. 

As a further comparison, we measured the radial abundance 
gradients in the youngest open clusters. The 24 open clusters 
with ages < 1.0 Gyr all lie within Rgc =13 kpc. The ra- 
dial abundance gradient for [Fe/H] is -0.05 ± 0.01 dex/kpc, 
a value in excellent agreement with the local Cepheids (-0.06 
± 0.00 dex/kpc). In the lower panel of Figure[32] we compare 
the radial abundance gradients, [X/H], for the young Cepheids 
and youngest open clusters (we only show data for those el- 
ements for which there were four or more clusters). The ex- 
cellent agreement in radial abundance gradients for the two 
sets of young objects validates our expectation that radial gra- 
dients represent a differential comparison within which the 
systematic errors cancel. 

One issue we can address with these data is whether the 
outer disk open clusters (Rgc > 13 kpc) share the same [X/H] 
abundance gradient as the inner disk open clusters (Rgc < 1 3 
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data, we superimpose the linear fits to the open cluster data. 



kpc). For open clusters with ages > 1.5 Gyr, we find that 
for all elements except La and Eu, the difference in gradients 
([X/H] vs. Rgc) between the local and distant samples is sig- 
nificant at the 3-cr level or higher. The same result holds when 
we consider open clusters with ages > 2.5 Gyr. Similarly, we 
can compare the abundance gradients ([X/H] vs. Rgc) for the 
local and distant Cepheids. In contrast to the open clusters, 
we find that for all elements, the abundance gradients agree 
within 1.8-<7. Furthermore, six of the 12 elements considered 
show agreement below the 0.9-er level between gradients for 
the local and distant Cepheids. In summary, the local and dis- 
tant Cepheids follow the same [X/H] vs. Rgc relation. There 
is, therefore, a fundamental difference in the behavior of ra- 
dial abundance gradients, [X/H] vs. Rgc, between the older 
open clusters and the Cepheids. 

Next, we compare the [X/Fe] radial abundance gradient be- 
tween the local and distant samples. For the open clusters with 
ages > 1.5 Gyr, the local and distant samples show the same 
[X/Fe] vs. Rgc gradient with the exceptions of Mn (2-cr), La 
(3-cr), and Eu (3-cr). When considering the open clusters with 
ages > 2.5 Gyr, with the exceptions of Mn (2-cr) and Eu (2- 
cr), the local and distant samples show the same abundance 
gradients. A similar comparison for the Cepheids reveals that 
with the exception of Al (2-cr), the local and distant samples 
show identical abundance gradients for [X/Fe] vs. Rgc at the 
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1.8-er level or lower. 

When we consider the trends between [X/Fe] vs. [Fe/H], 
our qualitative impression is that the local and distant open 
clusters follow the same relation. Similarly, the local and dis- 
tant Cepheids appear to follow the same [X/Fe] vs. [Fe/H] 
relation. Given the moderate ranges in [Fe/H] covered by the 
distant open cluster and Cepheid samples, these are qualita- 
tive rather than quantitative comparisons. For the linear fit to 
[X/Fe] vs. [Fe/H] across the full data range, we are willing to 
quantitatively compare the gradients between the open clus- 
ters and Cepheids. As before, we assume that the systematic 
errors in the analysis of a particular type of star will cancel 
such that we can reliably compare gradients between the two 
different classes of objects (red giants vs. Cepheids). We find 
that the open clusters and Cepheids share the same gradients 
(<3-cr) from the linear fit to [X/Fe] vs. [Fe/H] for all elements 
except O (3-cr), Mn (4-cr), and Ni (3-er). 

At the risk of over-interpreting these data, the results of 
this comparison ([X/Fe] vs. [Fe/H]) would suggest that the 
chemical enrichment histories of the older open clusters and 
Cepheids are surprisingly comparable. This similarity applies 
to the more metal-rich local samples as well as to the more 
metal-poor distant samples. Our conclusions would be un- 
changed had we increased our [Fe/H] ratios by 0. 1 dex and 
decreased our [X/Fe] ratios by 0. 1 dex for the open clusters in 
Paper I and in this study. Similarity in [X/Fe] vs. [Fe/H] re- 
lations between local and distant stars that span a large range 



in ages cannot be the result of the evolution of a closed sys- 
tem. It remains to be seen how such abundance ratios can be 
explained within a self-consistent chemical evolution model. 
One requirement may be the infall of gas onto the outer disk. 
However, that material must have experienced the same ratio 
of SNe II to SNe la, at a given metallicity, as the gas from 
which the local disk formed, both at early times (i.e., at the 
epoch of open cluster formation) and later times when the 
Cepheids formed. 

5. CONCLUSIONS AND FUTURE WORK 

Previous papers in this series studied the chemical compo- 
sitions of open clusters, field stars, and Cepheids in the outer 
disk. In this paper we present radial velocities and chemical 
abundances for a new sample of outer disk open clusters. This 
paper includes the first analysis of the outer disk clusters Be 
18 and PWM 4. We compiled a set of chemical abundances 
for a sample of 49 unique clusters drawn from our studies 
and from the literature. Using this sample, we studied trends 
between chemical abundances and distance, age, and metal- 
licity. 

We confirm the flattening of the metallicity gradient in the 
outer disk and that the outer disk open clusters are uniformly 
metal-poor with super-solar ratios for [a/Fe] and [Eu/Fe]. For 
some elements there are hints that the local [Rgc < 13 kpc) 
and distant (Rgc > 13 kpc) samples may have different radial 
[X/Fe] abundance gradients. 

We confirm that there are no significant trends between 
metallicity, or abundance ratios [X/Fe], and age (with the 
likely exception of the 5-process elements already noted by 
iD'Orazi et aT]|2009l and lMaiorca et alj|201 ll) . Compared to a 
sample of solar neighborhood field giant stars, we find that 
the open clusters share rather similar trends for [X/Fe] versus 
[Fe/H] for almost all elements. We quantify the linear trends 
between [X/Fe] and metallicity and find that the scatter about 
the mean relation, as low as 0.06 dex for Ni, is comparable to 
the measurement uncertainties for some elements. For other 
elements including Co, Zr, Ba, and La, the scatter about the 
linear trend is significantly higher than the measurement un- 
certainties which may suggest a real dispersion in abundance 
ratios. We note that for lines that are strong and/or affected by 
hyperfine structure (e.g., Ba), the measurement uncertainties 
may underestimate the true errors. Additionally, some ele- 
ments including Co and Zr are measured from small numbers 
of lines such that the measured dispersions are likely larger 
than the errors from stellar parameter uncertainties alone. Fi- 
nally, our analysis suggests that this inhomogeneous sample 
includes ^0.1 dex systematic offsets for some elements. 

The flattening of the metallicity gradient, differences in 
metallicity, and the enhancements in [a/Fe] (and perhaps 
[Eu/Fe]) suggest that the outer disk formed from gas with a 
different star formation history than the solar neighborhood. 
We reiterate that the individual a elements do not necessar- 
ily follow identical patterns. In particular, [O/Fe] and [Ti/Fe] 
are strongly enhanced in the outer disk relative to the inner 
disk while [Mg/Fe], [Si/Fe], and [Ca/Fe] show roughly con- 
stant ratios as a function of Galactocentric distance. Outer 
disk field red giant stars, which cover a more limited range 
in Galactocentric distance, share similar [X/Fe] and [Fe/H] 
vs. Rgc trends as the open clusters at the same range of dis- 
tances. When compared to Cepheids, the old open clusters 
(ages > 1.5 Gyr) show steeper [X/H] vs. Rgc trends suggest- 
ing that the Galactic disk grew via an "inside-out" process. 
The Cepheids and open clusters share very similar [X/Fe] vs. 
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[Fe/H] trends. Such similarity between samples of stars with 
very different ages cannot arise from the chemical evolution 
of a closed system. Understanding the chemical abundances 
of old distant open clusters and young distant Cepheids repre- 
sents a challenge for future chemical evolution models. 

Ultimately, a definitive statement about the origin and evo- 
lution of the outer disk requires a homogeneous analysis of 
larger samples of stars in larger numbers of clusters based on 
high quality spectra. Additional efforts should be undertaken 
to increase the samples of distant field stars and Cepheids with 
chemical abundance analyses for a more complete picture of 
the outer disk. Clearly, the comparison field stars also need to 
be analyzed on the same system as the open clusters. For ex- 
ample, the work by Przybilla, Nieva, and coll aborators (e.g., 
iPrzvbilla et a"T] 120081: INieva & Przvbillal 120 1 2|) on unevolved 
early type stars shows that very high precision abundances, 
10%, can be obtained for objects spanning a large range in 
distance. Extremely careful analyses may one day provide 
similar precision in abundances in open clusters to explore 



the origin and evolution of the outer disk. 
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APPENDIX 

A.l SIGNAL-TO-NOISE CONSIDERATIONS 

As noted in Section 4.1, the compilation of open cluster abundances includes data that span a range of S/N (and sample sizes 
per cluster). We now re-examine whether the trends between [Fe/H] vs. Rcc and [X/Fe] vs. Rcc change as we make cuts based 
on S/N. 

In Figure[33j we show [Fe/H] vs. Rcc again applying a linear fit to the local (Rgc < 13 kpc) and distant (Rgc > 13 kpc) samples. 
In this figure, we show (a) all data, (b) S/N > 50, (c) S/N > 70, and (d) S/N > 90. We find that the slopes and dispersions about 
the linear fit do not change substantially as we exclude data based on the S/N. One important point we highlight is that for the 
most stringent cut, S/N > 90, there is only one cluster beyond Rcc = 15 kpc. Panel (d) in this figure underscores the need to 
re-observe the most distant open clusters at higher S/N ratios in order to obtain more accurate and precise chemical abundance 
ratios. We also stress the importance of observing multiple stars per cluster, when possible, to increase the statistics and to ensure 
membership. Indeed, a similar figure could be generated taking into account the S/N, spectral resolution, and number of stars per 
cluster. 

In Figures [34]to [36] we plot [O/Fe], [a/Fe], and [Ni/Fe] vs. Rcc- For these representative elements, we again find that the slopes 
and dispersions about the linear fit to the local and distant samples do not significantly change as we restrict ourselves to only the 
highest S/N data. Nevertheless, we again argue that accurate and precise chemical abundance measurements are obtained from 
higher quality spectra and currently, there is a clear need to obtain such data for the most distant open clusters. 

A.l THE TRANSITION FROM INNER TO OUTER DISK 

iTwarog et al.l (11997b were the first to note that the radial abundance distribution, as traced by open clusters, shows a sharp 
transition near Rcc = 10 kpc. Numerous studies, including this series of papers, have sought to further understand the transition 
from the inner to outer disk. The location, and nature, of the transition can potentially provide crucial constraints upon the 
formation and evolution of the Galactic disk. While most studies have performed quantitative comparisons between the various 
(chemical) properties of the inner disk and outer disk open clusters, to our knowledge the transition region, or transition radius, 
is selected in most, if not all, studies through a qualitative visual inspection. In this appendix, we seek to quantify the location of 
the transition radius between the inner and outer disk. 

The data we shall use to quantify the transition radius are [Fe/H] vs. Rcc for the complete sample listed in TablefOl We assume 
that there are two regions, an inner disk and an outer disk. We also assume that each region can be described by a linear relation. 
Thus, there are five parameters to be determined: (1 and 2) the slope and intercept of the linear fit to the inner disk region, (3 and 
4) the slope and intercept of the linear fit to the outer disk region, and (5) the transit ion radius. 

We determine these parameters using the IDL MPFIT routine dMarkwa rdt 2009) which uses the Levenberg-Marquardt tech- 
nique for least squares minimization. One consideration is that an initial guess for the best parameters is required. On applying 
the MPFIT routines, this concern appears to be valid. We tested a number of initial guesses for the five parameters. The only 
parameter that showed a dependence on the initial guess was the transition radius. In the upper panel of Figure [37] we plot the 
initial guess (which ranges from 9 kpc to 16 kpc in steps of 0.25 kpc) and the optimum value. In this figure, there are two best 
solutions for the transition radius, 10.4 ±1.0 kpc and 15.3 ±1.3 kpc, and the optimum value depends on the choice of initial 
guess. (The Levenberg-Marquardt technique involves gradient descent. If there is a saddle point, it is therefore not surprising 
that the final solution depends on our initial guess.) We note that when using either of the two best solutions as the initial guess, 
their \ 2 values are almost identical, 49.86 (10.4 kpc) and 49.90 (15.3 kpc). In the upper two panels of Figure[38] we overplot the 
two solutions to our data. The first conclusion we draw from this quantitative analysis is that there does not appear to be a single 
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FIG. 33.— [Fe/H] vs. « G c for different S/N cuts. From top to bottom, the panels include (a) all clusters, (b) S/N > 50, (c) S/N > 70, and (d) S/N > 90. The 
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FIG. 34.— Same as Figure[33]but for [O/Fe]. 



value for the optimum transition radius if we do not place any constraints on the linear fits to the inner and outer regions (i.e., 
Case 1). Rather, there are two equally good solutions for the transition radius at 10.4 ±1.0 kpc or 15.3 ±1.3 kpc. 
The MPFIT routines include consideration of the errors for the Y axis ([Fe/H], for which we adopt a uniform value of 0. 15 dex), 
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FIG. 35. — Same as Figure l33l but for [a/Fe]. 
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FIG. 36.— Same as Figurel33lbut for [Ni/Fe]. 



but not for the X axis (Rgc)- In order to explore the effect of uncertainties in the distances, we adopt a Monte Carlo approach 
in which we replaced each distance with a random number drawn from a normal distribution of width 0.5 kpc centered at the 
Rqc of the given data point. We repeated this process for each data point in the sample and determined the optimum parameters. 
We repeated this process for 1,000 new random samples. As above, for each sample we trialed a range of initial guesses for 
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FIG. 37. — Optimum value for the transition radius as a function of the initial guess. Case 1 (upper) is when we do not require the two linear fits to intersect 
at the transition radius. Case 2 (lower) is when we require the two linear fits to intersect at the transition radius. The colors represent the \ 2 value for each data 
point. 



the transition radius (from 9 kpc to 16 kpc in steps of 0.25 kpc). In the upper panel of Figure[39] we plot the average optimum 
transition radius as a function of the initial guess. In this figure, the error bars represent the standard deviation of the distribution 
of 1,000 values. This figure again highlights that there are two preferred values which depend upon the initial guess. 

We note that thus far we have not required the two linear fits to intersect at the transition radius. We now consider how the 
results change if we constrain the two linear functions to intersect at the transition radius (i.e., Case 2), and this may be a more 
appealing way to describe the behavior of metallicity with Galactocentric distance. In this scenario, there are only four free 
parameters (the slope and intercept of the fit to the inner region, the transition radius, and the slope of the fit to the outer region). 
Once again, we test to see whether the optimum value for the transition radius shows a dependence on the initial guess. In the 
lower panel of Figure [37] there is a suggestion that there are again two solutions which depend on the initial guess, 11.4 ± 0.9 
kpc and 12.5 ±1.1 kpc. In the lower panels of Figure[38] we overplot these two solutions and note again that the \ 2 values (when 
adopting the optimum values as the initial guesses) for the two solutions are very similar, 50.48 (1 1.4 kpc) and 51.53 (12.5 kpc). 

Finally, we apply the Monte Carlo approach described above to Case 2. In the lower panel of Figure[39]we plot the optimum 
transition radius as a function of the initial guess where the error bars represent the standard deviation of the distribution of 1,000 
values per initial guess and the color-bar represents the average \ 2 - In contrast to Case 1, the Monte Carlo simulations suggest 
that for Case 2 there is a single best solution for the transition radius of 12.1 ± 0.7 kpc. The second conclusion we draw from 
this analysis is that if we require the linear fits to the inner and outer regions to intersect at the transition radius, there is a single 
optimum value for the transition radius of 12. 1 ± 0.7 kpc. 

Finally, we note that within the uncertainties, the slopes and linear fits we obtain from the MPFIT routines for the inner and 
outer regions agree with the slopes we measure in the main text based on our qualitatively selected 13 kpc transition radius. 
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FIG. 39. — Same as Figure |37l but based on a Monte Carlo simulation for N = 1,000 realizations. In each realization, the Rqq for a given data point was replaced 
by a random number drawn from a normal distribution of width 0.5 kpc centered at the original value. The data and error bars represent the average value and 
dispersion from the 1 ,000 realizations respectively. The colors represent the average \~ at eac h data point. 
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